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In  the  special  problems  considered  in  sections  17  to  21,  the 
dependence  upon  the  time  has  been  precipitated  out,  either  by 
assuming  the  motion  steady  in  a  moving  coordinate  system  or  by 
assuming  a  harmonic  dependence  upon  the  time.  In  this  section  we 
shall  consider  motions  in  which  the  displacement  and  velocity  of 
the  surface. are  specified  at  some  instant  of  time,  say  t  -  0,  the 
motions  of  any  solid  boundaries  are  given  for  each  instant  t  >  0 
(exce.pt  at  the  very  end  where  freely  floating  bodies  are  considered) 
and  the  pressure  distribution  over  the  free  surface  la  a  given  func¬ 
tion  for  t  >  0. 

It  is  not  usually  possible  in  the  most  general  situations  to 
give  explicit  solutions  for  such  problems.  However,  Volterra  [1034] 
has  proved  a  uniqueness  theorem  and  has  shown  how  to  reduce  the 
problem  to  finding  an  appropriate  Green's  function.  His  results 
were  later  rediscovered  and  extended  to  a  wider  class  of  problems 
by  Finkelstein  [1957]  (see  also  Stoker  [1957,  ch.  6]).  However, 
the  use  of  Green's  functions  for  initial-value  problems  extends  back 
even  earlier,  at  least  to  the  papers  of  Hadamard  [1010,  1016]  and 
Bouligand  [1913].  These  theorems  are  discussed  in  section  22m. 

One  of  the  classical  problems  in  this  category  is  associated 
with  the  names  Cauchy  [1827]  and  Poisson  [1815].  In  this  problem 
the  fluid  is  infinite  in  horizontal  extent,  without  obstructions, 
and  either  infinitely  deep  or  of  uniform  depth  h.  At  the  initial 
instant  t  -  0,  the  form  of  the  surface  and  its  vertical  velocity 
are  given  and  one  seeks  the  subsequent  motion.  Such  problems  have 
already  been  discussed  at  some  length  in  section  15.  However,  in 


3  If 


that  section  interest  centered  upon  investigation  of  certain  aspects 
of  the  subsequent  motion  rather  than  upon  obtaining  the  solution. 

In  addition,  the  treatment  of  that  section  was  limited  to  two-dimen¬ 
sional  motion,  although  the  methods  could  have  been  extended  to 
three-dimensional  motion. 

The  history  of  this  problem,  including  an  exposition  of  the 
methods  used  by  various  authors,  is  Included  in  a  paper  by  Risser 
[1924,  pp.  113-144].  Another  expository  account  can  be  found  in 
Vergne  [1928,  ch.  I].  The  problem  is  discussed  here  in  section  22 

In  section  22  ft  several  special  initial-value  problems  are 
discussed. 


22  cur , 


Some  gen< 


neorems 


Let  the  fluid  be  bounded  by  the  free  surface  F,  fixed  surfaces 
S  and  the  surfaces  of  a  finite  number  of  bodies  of  bounded  extent 
undergoing  specified  mottoes  of  small  amplitude  about  equilibrium 
positions  Sm  (see  Figure  28).  Let  the  pressure  dlatrlbution 
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on  the  fro©  aurfac®  F,  also  be  a  given  function  p(x,a,t).  Further¬ 
more,  at  time  t  »  0  let  the  Initial  displacement  and  vertical  velo¬ 
city  of  the  free  surface  be  given  functions 

7,'*, 5.o)(  Vt  -  {*.},«>;.  .  (M-1) 

The  boundary  conditions  to  be  satisfied  by  the  velocity  potential 
(x,  y ,  z,  t)  are  (see  11.1) 

4*n  *  ( i)  on  Sm/  (22.2) 

ft  (*.0.  ],o)  *  -  <j  *1  (*'  l>°)-  ]'  T* 5'  ^  ° n  F/ 

^  (x,0,  ),o)  =  *?f  )'°)  on  F  . 

Here  F  means  that  part  of  the  plane  y  ■>  0  occupied  by  fluid  when 
everything  is  at  rost.  In  aduitlon,  it  will  be  assumed  that,  for 
each  t,  there  Is  a  bound  B  and  a  distance  rQ  such  that  |^|  ,  |  ft\, 

(grad  f\  and  |grad  ft\  are  each  less  than  B  for  x2+y2+ as2  ?  rQ2. 

Let  us  now  suppose  that  it  is  possible  to  find  a  source  func¬ 
tion  G  of  the  following  nature: 

1,  5  i  i.  v)  -  ±  ,  S,  1,  i ;  I,r), 

/nn  *%  \ 

\ w / 

where  as  usual  r2  -  (x-f|)2+(y-  £)2+  '.a-  p2  and  H  is  harmonic  for 
y  ^  0,  and  in  addition  G  satisfies 


6  (*>  1.)  ;  Jf  ,  *1,  5  ;  t  )  -  (*•  1')  i  ?  *  5  ;  t,  t)  *  <?, 


(22.4) 


<ji  j  5,  7,  ^  ;  t,  r)  «  o  -f°r  te'I'V  onS  anaL 

{or  all  i 


This  function  has  already  been  constructed  in  two  cases.  If  there 
are  no  fixed  boundaries  and  the  fluid  is  lm  Initely  deep,  the  func¬ 
tion  defined  in  (13.49)  satisfies  the  conditions  after  slight  modi¬ 
fications  .  replace  (a,b(c)  by  (£,  ^  ,  *  ) ,  set  m(t)  -  1,  and  extend 
the  definition  of  (j>  (of  13.49)  to  negative  t  by  y,z,-t)  - 
^(x,y,SB,t).  Then  we  take 

CU,<j,3  >  ^,7.  ?  i  t,  r)  ~  tuy] ;  t-r)  ■ 

Similarly,  the  function  defined  in  (13.53)  allows  one  to  construct 
G  when  the  fixed  boundary  consists  of  a  horizontal  bottom  at  y  •  -h. 
For  the  first  G.  Finkelatein  (1957,  Appendix]  has  shown  that  G  is 
0(R"2)  and  Gr  and  Gy  are  0(R~3)  as  R*-o,  where  R2-(x- J  )2+<z-  J)2; 
for  the  second  G.  Finkelstein  [1957,  fa]  has  shown  that  0,  GR  and 
Gy  are  ♦  €)%)  *s  R  o©  for  arbitrary  t  >0  , 

Now  apply  Green's  theorem  to  the  functions  <f>t  and  G  and  the 
region  of  fluid  bounded  by  the  surfaces  Sro,  the  fixed  boundaries 
S  and  a  large  sphere  Jl  of  radius  f  and  center  at  the  origin, 
where  f  is  chosen  large  enough  to  include  all  the  surfaces  S,q. 

Only  parts  of  F,  S  and  J2  will  serve  as  bounding  surfaces,  and  we 
shall  call  these  parts  F',  S'  and  J2. ' ,  respectively.  Then 


fvs 


7(22.5) 


where  v  is  the  exterior  normal.  The  right-hand  side  is  actually 
independent  of  T  since  T  enters  only  through  the  function  H  which 
is  harmonic.  The  integral  over  S'  vanishes  since  both  </>  n  and  GQ 
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arc  aero  on  S.  We  shall  assume  that  th«  behavior  of  G  and  <f  a# 

1  ♦  *o  la  such  that  ths  Integra l  over  SL  '  vanishes  as  f-aao1.  If 
the  fluid  is  of  bounded  extent,  the  situation  considered  by  Voltsrra 
[1934],  this  presents ,  of  course,  no  difficulty.  In  ths  two  cases 
for  which  G  has  been  given  above,  it  has  been  shown  by  Finite  1st® In 
that  this  is  true.  For  finite  depth  the  proof  presents  no  diffi¬ 
culty  one®  ths  estimates  for  G  ars  obtained;  for  infinite  depth  the 
analysis  is  mors  troublesome  and  ws  ra fsr  to  his  paper  or  to  Stoker 
[1057,  pp.  103-4]  for  proof.  After  letting  f  -*■*>»  one  then  hae 

Ij*')']'*)  m4T  ^7  ~  S  ^1***4$  f(fr  Av  ^d6- 

F  1  1  v  (22.6) 

■*  1 

In  the  integral  over  F  we  may  replace  by  -g  Gtt  btcau*®  °* 

the  boundary  condition  at  F.  Now  interchange  t  and  V  and  integrate 
with  respect  to  T  between  limits  0  and  t.  This  gives,  following 
an  integration  by  parts, 

■to  n  n  -  f  «>5  ■ ■  °>  |j/fc  *„  *  $  it  j|‘-  fk  i>7 *■ J?  ^ 

F  0 

<M-T) 

-ii  //  {c  a, °,i, * .  i ,  j  i  * .  *. )i<j  t\. . °>  i-  *)*  \  h' •  > e)  Gt r  "j  *• 

F 

-  G("  v'Q  ‘,0)~  ^  o)Gt(,..;  0t  i) 

+  ^  (■  ■  •;  T>  0  Pi  (tt  r)  dT  |  d  d ^  ^ 

where  I  stands  for  the  last  integral.  (Gt  always  represents  the 
derivative  with  respect  to  the  seventh  variable).  Recalling  the 

. 'Hi NiMi mm . me tMumma  Ju •  M"«im •  '  “iU* * •  l"lllllM"' w 1  ■ 
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properties  of  G  In  (22.4),  one  finds 


t 

(22.8) 


?lx‘1  •)'*)  -  ?(*•'},}, °)+fr  (J  f-C(flo,5l*.l'}i°.t)ll(f,l,o) 

*  **  1‘‘ *■  i‘  a ; *  Ofo  (h  W+h  r(ll,°)]*fa  fa  (l,°,  ?,•*,,  .},r,t)  ■ 

ri(t,l,T)dT}*l  di+- 


-  o)  *■&  fsj-cdo.i,  O,  t)  ^  ($ ’,14+  Gt  (f.0,1,  ot)<j(l  j,o) 

.  ^  faao.i,  jc.^j^.t'rfrsxpTjd^i  +  it 

where  (x, y,z,o)  is  determined  up  to  an  additive  constant  aa  the 
solution  to  a  Neumann  problem,  since  -n(x,y,z,o)  is  given  on  all 
boundaries  and  bounded  at  infinity.  In  the  integral  I  we  note  that 


^  tn  "  vn '  (t)  is  known  on  Sg, f  but  ^  t  is  not. 

If  there  are  no  moving  bodies  In  the  fluid,  then  the  integral 
I  Is  not  present  and  la  determined  by  the  initial  displacement 
and  velocity  of  the  free  surface  and  the  given  pressure  distribution 
over  it.  This  is  Volterra’s  result  as  extended  to  unbounded  fluids 
by  Finkelsteln.  If  surfaces  Sm  are  present,  one  may  still  use 
(22.8)  to  derive  an  integral  equation  in  the  same  way  that  (16.13) 
was  derived.  For  as  (x,y,z)  is  made  to  approach  a  point  (x0,y0,z0) 
of  Sm, 


Sm 

+4TT  ■*’>  h  da 

lef,  Kellogg,  Foundations  of  potential  theory,  Springer,  Berlin, 
1929,  p.  167).  Thus,  after  carrying  out  the  integration  with 


. . 


i 

I 


respect  to  T  ,  one  has  an  Integra!  elation  tor  f  tor 

each  value  of  t  7  0  which  say  be  ueed  to  find  the  value  of  ^ ,  and 
hence  on  the  surface  Sn,  providing  that  the  integral  equation 

can  be  solved.  One  nay  then  use  (22. 8)  to  determine  d  (x,y,x,t) 
for  ail  values  of  (x,y,z)  \  the  x  i.  The  integral  equation  has 

the  same  appearance  as  (22.8)  except  that  the  first  two  te  rata  have 
coefficients  §  and  (x,y,z)  is  understood  to  be  a  point  of  Sm,  This 
further  extension  of  Volterra’s  analysis  is  also  due  to  Finksisteln. 
Uniqueness  of  (x,y,z,t),  at  least  up  to  an  additive  constant, 

I  ' 

may  be  proved  as  follows.  Let  ^>\  and  <p  j  be  two  solutions  satisfy¬ 
ing  the  boundary  conditions.  Then  <P  •  <P\  **  $>  2  aatlsfles  (22.8) 
with  f,F,  p  and  Vra  all  Identically  sero,  l.e. 

f  -h(dT 

!  If  we  assume  that  Gn  is  0(R_1“f)  as  R  -►  ®o ,  then  <f>  t  and  grad  <f» 

will  have  the  same  behavior  and  the  Integrals  we  shall  write  below 
may  be  shown  to  exist.  As  has  been  mentioned  above,  Gn  vanishes 
much  quicker  than  is  required  in  the  cases  when  the  fluid  is  lnfi- 

1 

nitely  deep  and  when  the  fixed  surface  consists  of  a  horizontal 

bottom;  if  the  fluid  is  bounded  in  extent,  no  such  condition  is 

I 

necessary  to  make  the  integrals  converge.  < 

Consider  then,  following  Volterra, 

11  “  i  T{  H  ^  eL6 

F 


=  -  //  <h  tr,  d  e 

F  f  S  + 


I 
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wince  $  R  vinlihes  on  S  and  SB.  Now  apply  Green's  theorem  and  de¬ 
note  the  volume  occupied  by  fluid  by  T: 

J2  -  -  Jf(  i  dr  -  -  i  2-t  ff farad  4>l)dT 

r  t 

Hence 

ii  jffj  **  *  Tj  -  ° 

and 

^  da  +  JJfafrlKd.  »  const  , 

(22.0) 

Since  <f>  n  -  0  on  F,  S  and  SB  for  t  -  0,  ^(x,y,z,o)  -  C,  a  constant; 
hence  grad  ^  -  0  for  -  -  0.  Also  ^t(x,y,z,o)  «  0.  Hence  the  con¬ 
stant  in  (22.9)  is  zero  and  i  t  and  grad  ^  vanish  for  all  t.  Thus 
(x,y,s,t)  -  const,  and  the  solution  of  the  initial-value  problem 
la  dstsrminsd  up  to  a  constant. 

22  ^  .  The^auch^^ 

In  this  classical  problem  of  water-wave  theory,  the  pressure 
over  the  free  surface  is  constant,  say  zero,  the  fluid  is  Infinitely 
deep  or  bounded  below  by  a  horizontal  bottom,  no  obstructions  are 
present  and  the  initial  displacement  and  velocity  of  the  free  sur¬ 
face  are  given.  The  two-  and  three-dimensional  cases  will  be  se¬ 
parated  in  order  to  illustrate  different  methods  of  approach. 

IIl£g£..dl.n»pg i ons .  The  velocity  potential  may  be  obtained 
directly  from  (22.8)  after  setting  jo  (x,z,t)  and  I  equal  to  zero. 
However,  the  explicit  expressions  for  G  and  Gt  are  needed.  As  was 

ll 
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noted  Is  unction  22 Athene  can  be  written  down  immed lately  from 
(13.49)  for  infinite  depth  and  (13.53)  for  depth  h.  The  reiiul, 
expressions,  after  aattlag  7  -  0,  art  as  follow*: 
infinite  depth: 

i.o.i  zfc‘-c~(qv  t)]*^  J.ck  K)dk, 


(22.10) 


Qt  (*,}}  i  f,  0 ,  $  ;  o,t )  -  -  2[  Su*.  ( iTyT  t)c  ?  J„  fa  R)  d  k  j 


depth  h: 


(22.11) 


(*, yyt f,o, j f 0, t) -  -2 nx  (J^TZIU t) 


where  U  1 «  ( je  -  f,  )V  t  (]  *-  S)  . 

There  etill  remains  to  find  ^<x,y,z,o)  where 


4>)(*,<>t},0)  •  ^4  (*,  5,0) 


and 


^  ^  ( »x ,  ^ j , »)  *  o  °r  ^  ^  (x  r  j  *  °)  9  0‘ 

(O  ■*  —  w£» 

The  solution  of  these  two  problems  is  well  known: 
infinite  depth: 

tO  ,  c  ,  .  *°  T ku 

a,  t  M  x  «i.i  if  -— — -~rrpr  ab  u?  e 

(22. 12) 
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depth  h 


•*  -■ 

J-(<e 


(22.13) 


where  E  is  defined  as  above. 

Subatitutinf  the  several  expressions  In  (22.8),  one  obtains 
the  expressions  for  the  velocity  potential: 
infinite  depth: 
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iw  J[l(tllo)dld^6^  &AMS  6  £  J«  (U  £ -)cL 
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depth  h: 


oO  r>0 

lo)d?d}[  6»*tTl(±*)4k 

**»  ffol J!ffWu  ^  J'f4 
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The  equations  describing  the  free  surface  are  as  follows : 
infinite  depth 
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depth  h: 
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It  has  been  shown  by  (Cochin  (1935]  thst  tbs  Integrals  with 
rsapsot  to  k  in  (82.16)  can  bs  evaluated.  Consider  the  integral 

00 

K  *  /  ®“  .U**,  or  f  J*  ^4  R  ^  <a(  4  /  , 

*  * 

(22. 18) 


Then  the  first  Integral  with  respect  to  k  in  (22.16)  is  ~Itt  and 
the  second  one  Is  To  evaluate  K  make  first  the  following 
change  of  variables: 


Mn  L%  jMk  mu 

i  non 


Jt  -  fc  R  ,  <j  t  U  R 


k  -  j  (tf1) 
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(22. 19) 
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In  the  first  integral  let  u 


v<*  +  W  ,  in  the  second  let  u  -  vM-iq. 
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Then 
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•nd,  after  acttlng  u  -fT  uJ 

K'^u[  <T.(*  £‘*>9) 

Finally ,  from  an  identity  In  Watson 'a  Bessel  function*  [|5„43,  eq 
(1)]|  one  finds 


fC( u>; «  & 


(22.20) 


In  order  to  use  the  result  in  (22.16)  one  needs  the  first  three 
derivatives  with  respect  to  t.  Since 

"s  .  rr  3 

51  ■  *»t  *  ' 

|he  derivatives  can  be  computed  by  taking  derivatives  with  respect 
to  to  and  multiplying  by  an  appropriate  factor.  After  some  rather 
tedious  computation  one  finds 

£  K<">  -  -  4  >'£[<]<  a  »)Jt  a  ~'j- tyt  »‘Mt  a  «vj, 
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K(w)  -  -  /  J2  (  i  x},(i^)-  «//'  (  i  '+>')  JlU  *x) 

v*  7*  tL,  *‘  4  '  **%  *•#  ^ 

"  \J  ,  ( {  u l)  tl  ^  (  t >Ai  )  J  J 

These  are  Kochin 's  formulas,  but  derived  soma  what  differantly  from 
hia  original  paper;  still  another  derivation  nay  ba  found  in  Kochin, 
libel  and  Roza  (1048,  ch.  8,/  21),  Similar  formulae  for  (33.17) 
do  not  seam  to  have  baan  dlacovered. 

It  should  ba  noted  that  in  the  find  form  of  (33. 16)  the  depen¬ 
dence  upon  t  la  through  the  dlmena  ion  lean  variable  to  a  •  gta/4R, 
Hence,  if  one  examinee  the  contribution  to  the  surface  profile  from 
a  given  locality,  aay  the  neighborhood  of  (J  ,  £),  then  a  given 
phase  of  this  contribution,  say  a  maximum,  will  be  described  by 
gta/4R  -  const.  ;  i . e . ,  the  phaee  is  moving  away  from  ( £  ,  £ )  with 
constant  acceleration  proportional  to  g.  The  amplitude  of  the  con¬ 
tribution  is  modulated  by  either  r_3/2  or  g-2  according  as  one  is 
considering  the  first  or  second  summand  in  (22.16).  Kochin 's  1039 
paper  is  of  some  methodological  Interest  inasmuch  as  he  started  his 
analysis  with  dimensional  considerations.  This  method  will  be  in¬ 
troduced  for  the  two-dimensional  case. 

One  may  obtain  without  great  difficulty  series  expansions  for 
the  k-lntegrals  in  (32.14)  and  (22.16),  as  was  first  done  by  Cauchy 
and  Poisson,  le  refer  to  Lamb1*  Hydrodynamics  [1032,  §  255]  for  the 
derivation  and  exact  expressions.  They  can  also  be  derived  from 
the  the  known  expansions  for  Jj,  etc.,  as  can  asymptotic  expressions 
for  large  to  .  One  may  also  carry  out  an  analysis  of  the  changing 
shape  of  the  surface  profile  following  the  methods  of  section  15. 

It  is  evident  that  one  can  solve  explicitly  other  similar 
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initial-value  problems  for  which  the  Green's  function  can  be  given. 

For  example,  the  method  of  Images  allows  one  to  give  an  explicit 
solution  for  various  cases  when  vertical  walls  are  present  as  boun¬ 
daries.  Such  cases  have  been  considered  by  Elsser  [19291.  The 
Cauchy-Poisson  problem  In  the  presence  of  a  vertical  half-plane, 
a  -  0,  17O,  has  been  treated  by  Boiko  [1938],  but  by  more  complex 
methods. 

Two  dimensions.  Rather  than  repeat  the  methods  used  for  three- 
dimensional  motion,  we  shall  introduce  a  method  making  use  of  the 
complex  potential  and  thus  special  to  two-dimensional  motion.  It 
Is  analogous  to  ths  method  used  in  deriving  (13.28). 

Let  {($ it)  •  be  the  complex  velocity  potential. 

The  initial  conditions  will,  be  taken  in  the  form 

“J  0-0“  ?<*•)/  •  (22.22) 

Let  us  consider  infinite  depth  first.  For  t  ?0  we  assume  that 
f (z, t)  is  regular  and  if )  <  M(+),  /  <  tl(p)  for  y  <■  0  and  that 

both  f  and  f  ft  approach  zero  as  y  -*•  -  ao.  Consider  now  the  func¬ 
tion 

fttO'O  +  1  j  /  (}fi)  ■ 

(22.23) 

From  the  assumptions  about  f  it  follows  that,  for  t  ?0,  G(z,t)  Is 

that  I  Cl  <  -6^ for  ^  <  °  and  that  G-*0  as 
«*■  -oo  .  Moreover,  it  follows  from  the  condition  at  the  free 
surface,  (11.5),  that  Re  G(x-  C o, t)  -  0.  Hence,  the  definition  of 
G  may  be  extended  into  the  upper  half-plane  by  defining 
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But  then  since  G  Is  regular  and  bounded  in.  the  whole  finite  a-plane, 
it  follows  from  Liouvi lie’s  theorem  that  G  -  const.,  the  constant 
must  equal  aero  from  the  assumed  behavior  as  y-r-  **> .  Hence  the 
fundamental  differential  equation  fur  the  Cauchy-Poisson  problem  in 


two  dimensions  is 


f  ;« 


(22. 25) 


an  observation  usually  credited  to  Levi-Clvita  [of,  Tonolo,  1913]. 

Let  us  now  find  the  analogous  equation  for  finite  depth.  The 
function  f  will  be  assumed  regular  in  the  strip  0  ?  y  7  -h.  The 
boundary  condition  on  the  bottom  is 


^  T  ix' 1  t)  ° 


(22.26) 


Hence  f  may  be  extended  analytically  into  the  strip  -2h  < y <  -h  by 


defining 


■f  (x  -  t  (  b  ^  Q)  ■  -f  ( *-  *■  1  ”j  )  /  0  *?  ^  ^ 


(22.27) 


We  may  also,  as  before,  extend  the  function  ftt*iKf 1  lnt0  the  s,triP 
h^y^O.  The  condition  Re  [  ftt+18J  j  “  0  for  Y  m  0 
Re  ]ftt~igf'}  -  0  for  y  «  -2h.  Hence  the  function  ftt-iBr'  can  b® 
extended  by  reflection  into  the  strip  -3h ^  y £ ~2h.  Now  consider  the 


function 
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(22.28) 

*  j  -fit  ( ]  *  ‘  ^ / 1)  + 4 j|  {  /  J *  ‘^r  ^  j +  |Vh  f  J  "  t  ^  j  ”  ^  0  1  ^  ^  J . 

A#  a  result  of  the  various  extended  regions  of  definition,  one  nay 
verify  easily  that  H  is  defined  for  all  *  in  the  strip  ~2h<y  <0 
and  is  regular  there.  Moreover,  it  follows  that 

H  (X~  i  k,  t)  •<?  )  (22.20) 

for  from ©2.27)  it  follows  that  the  two  pairs  of  sumn.ands  in  the  first 
form  of  (22.28)  are  real,  for  z  -  x-ih,  whereas  from  the  boundary 
conditions  at  y  **  0  and  y  **  -2h  it  follows  that  the  terms  in  curly 
brackets  In  tne  second  fora;  of  (22,28)  have  zero  real  partis.  Since 
H ( z , d )  is  regular  in  the  strip  0  7  y  t  -2h  and  vanishes  on  y  -  ~h,  it 
Must  vanish  identically  in  the  strip.  Hence  we  have  the  following 
differential-difference  equation  of  Cisotti  (1018]: 

Wj  +it m °-  (22*30) 

Let  us  now  turn  to  the  solution  of  (22.25)  with  initial  condi¬ 
tions  (22.22).  We  shall  follow  closely  an  exposition  of  Sedov’s 
[1048,  1957].  However,  the  idea  of  the  derivation  is  Kochln’s  [1035] 
and,  in  fact,  really  goes  back  to  Tonolo  [1913],  The  use  of  dimen¬ 
sional  analysis  can  be  extended  to  the  three-dimensional  problem; 
this  was  also  done  by  Kochio. 

We  first  remark  that  the  initial-value  problem  can  be  solved 
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by  solving  it  for  two  special  cases  of  (22. 22) ,  namely,  first  with 
n  ( x » 0 )  •  0  ana  then  with  (*,«>>  *  <>•  Tbs  su m  of  these  two. 

solutions  will  satisfy  (22.22).  Neat  we  note  that  *j  ( a , 0)  has  the 
dimension  "length"'  and  ^t(x,0)  the  dimension  "velocity",  and  that 
the  solution  f  in  each  of  the  t*  >  initial -value  problems  will  be 
proportional  to  some  typical  parameters  associated  with  *7  U,0)  or 
{ x  1 0 )  1  respectively.  Let  us  suppose  that  a  is  such  a  parameter 
with  dimension  Lf"  T<1  and  that  f  is  proportional  to  a.  Sine©  f  has 
dimension  L2T-1  and  g  has  dimension  LT~2,  the  Pi  theorem  of  dimen¬ 
sional  analysis  then  states  that  f  can  be  expressed  as  follows: 

Sl  <*■}*  ^  *  (  ~fr  )  t  (22.31 


where 
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(22.32) 


(The  factor  i/4  in  the  argument  of  is  chosen  for  later  conven 
ience . )  Now  substitute  (22.31)  into  (22.25).  One  finds  after 
iome  computation  that 


*f+t  + 


)  (22.33) 


where  T**  i  ^  M  J} 
ting  the  expression 


.  The  differential  equation  obtained  by  set- 
in  square  brackets  equal  to  aero  determines 


^ In  terms  of  confluent  hypergeometric  functions: 


■t(5)=AtF,(-*.ii3)+-&  VfatCi-*, \>S) 


(22.34) 
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From  this  it  follows  that 
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Remembering  that 

,pt  (v,  r jo)  -  *  ,  fid,  fi°)*  Vf, 

one  may  easily  derive  the  following' 

-  A  i 

,  ,  fi  «x~l  (22.36) 

t  M'°)  *  A  f,  ()<°) w  A  a  °c  ^  }  , 

.  I  ^  (*  -  * 

f*(V°)  =  &  *  i 5  ^  2 

The  solution  (22.33)  may  be  further  generalized  by  replacing  t  by 
t~t0  and  z  by  z-xQ  (l.e.,  by  a  different  choice  of  the  dimensionless 
variable  f  ).  One  may  then  further  superimpose  these  solutions. 

For  the  purpose  at  hand  it  will  be  sufficient  to  retain  t0  "  0. 

Then  we  may  font  the  solution 

nnO 
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One  finds  from  (22.30)  that 
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• 1  J»s, 


Let  us  now  make  some  special  choices  of  a,  and  hence  of  ec  and 
/£.  As  a  parameter  describing  the  initial  profile  of  the  surface 
take 


a 2  -  j  1  (*>o)  dx  ■ 


(22.39) 


as  a  parameter  describing  the  initial  velocity  distribution 
take 

iwO  X.  ^ 

a,  =  f  dx  (  y  (l,o)*l  . 
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(22.40) 


Then  »i  has  the  dimension  l‘lT . *,  corresponding  to  Q(|  -  -l,  j  ft"  0, 

and  a 2  the  dimension  L2,  corresponding  to  o  2  “  /-^  *"  I*  With 

these  choices  of  (X  ^  and  (X  2  in  (22.37)  we  take 
_  .  ** 
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Then  the  last  two  equations  of  (22.38)  become  (after  an  integration 
by  parts  in  the  first  one) 
h  .  1 
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From  the  Plemel J-Sokhotskii  theorem  we  have 

*0 

/  (X  -  <1 0,  oj  »  -  t  F(x;  +  -  d  x, 

—  oO  0 


(22.42) 


Thus  the  Initial  conditions  {22.22}  are  satisfied. 


There  remains  to  paint  out  that  for  the  special  choices,  of 
(X  |  »  -1  and  «2  -  -  J  the  corresponding  confluent  hyper  geometric 
functions  In  (22.35)  may  be  expressed  in  terms  of  Fresnel  integrals 


or  integrals  of  these.  In  fact,  if  we  write  (22.37)  in  the  form 
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then 


(22.44) 


On©  should  also  consult  the  discussion  in  Lamb’s  Hydrodynamics 
(1932,  £  238,  239],  where  graphs  are  given  which  display  the  behavior 
of  the  surface  profile  corresponding  to  an  initial  elevation  con¬ 
centrated  in  the  neighborhood  of  one  point,  i.e.,  essentially 

~  1  t 

-g  i2  2t(*-i0#  t)  *  and  to  a  concentrated  impulse,  i.e.,  essentially 
-g~l  It  it(x-tO,  %) ,  However,  general  aspects  of  the  development  of 
the  surface  profile  have  already  been  discussed  in  section  15  oc. 

It  should  be  noted  that  the  velocity  potential  (22,37)  repre¬ 
sents  a  much  wider  class  of  time -dependent  gravity-wave  motions 
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than  does  (22.43).  The  initial-value  problems  corresponding  to 
other  valves  of  aC  have  been  determined  by  Sedov  [1948]  but  the  dis¬ 
cussion  will  not  toe  repeated  here. 

A  class  of  solutions  of  (22.30)  analogous  to  that  found  by 
SedO'V  for  (22.25)  does  not  seem  to  have  been  given  in  the  published 
literature.  Cisottl  j 1920]  expands  l(z,t)  in  a  power  series  in  t» 
thus  replacing  (22,30)  by  a  recursive  set  of  difference  equations.. 
We  refer  to  the  original  paper  for  his  discussion  of  this  set  of 
equations. 

In  (15.22)  we  have  already  given  the  velocity  potential  ard 
surface  profile  corresponding  to  a  given  initial  profile;  the 
derivation  wan  based  upon  a  Fourier  analysis  of  the  initial  profile 
and  the  result  was  valid  for  either  finite  or  infinite  depth.  The 
same  procedure  may  be  used  for  an  initial  velocity  distribution. 

The  combined  result  for  the  complex  velocity  potential  and  surface 
profile  is  given  below  in  such  a  way  as  to  include  the  possible 
presence  of  surface  tension: 
infinite  depth; 

-f ( J  i i )  -  ~  j  d k  j d  |  l  j’  <5 (  k)t}{ 1 0)  w.  at  t^{(  6  i  j  t  ^ ^  ^ 
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where  &  ,c  |fl';J^E  ,f  ^  /  J-O'vX  ^ 

If  t  -  0,  the  coefficients  of  «j>t(  f  ,0)  In  the  formulas  for^(x.t) 
reduce  to  <5  "l.  then  T  -  0,  (22.45)  Is,  of  course,  another  form  of 

(22.43). 

It  has  already  been  Indicated  in  section  15  that  the  Cauchy- 
Poisson  problem  can  also  be  solved  for  two-fluid  problems. 

Sretenskil  (10551  has  investigated  a  further  generalization  in  which 
the  two  fluids  are  each  flowing  with,  constant  velocities  for  t<0 
and  then  when  t  -  0  a  disturbance  is  suddenly  created  at  the  inter- 


22  jf  . 

It  is  possible  to  solve  a  twn.be r  of  initial-value  problem* 
either  by  using  equation  (22,8)  or  by  using  the  time -dependent 
Green  functions  (13.49)  or  (13.53)  directly.  The  special  situa¬ 
tions  treated,  below  fall  roughly  into  the  following  four  cate¬ 
gories:  wave  motions  resulting  from  a  pressure  distribution  sudden¬ 
ly  brought  into  existence  at  time  t  -  0,  waves  resulting  from  a 
body  set  into  motion  at  time  t  -  0;  waves  resulting  from  an  under¬ 
water  explosion  or  a  sudden  movement  of  the  bottom  (tsunamis);  and 
waves  resulting  from  an  initially  displaced  freely  floating  body. 

Time-de pendent  pressure  distributions.  Suppose  that  the  fluid 
is  undisturbed  for  t  <  9  and  that  starting  with  t  -  O'  the  pressure 
over  the  free  surface  is  a  given  function  p  The  consequent 

motion  of  the  fluid  may  be  easily  obtained,  for  this  is  just  the 
problem  formulated  in  (22.2)  if  we  put  )”  Y/JX']t'°)*0’  Formu~ 

la  (22.8)  then  gives  the  velocity  potential  in  the  form 


*  F  0 

In  the  two  situations  for  which  explicit  Green’s  functions 
have  boen  given,  equations  (22.10)  and  (22.11),  we  may  give  expli¬ 
cit  solutions  for  ^  : 

Infinite  depth; 
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where,  as  usual,  R  •{■X-  Jf  )  f  -.J y  *. 

The  velocity  potential  for  a  moving  pressure  distribution  is 
obtained  from  these  expressions  simply  by  letting 

p(f,},T)=  f.d-cr,  It) 

It  l  . CT,  Z/T)-p.(f-cT,pc*S<FZ  the  resulting  f  is  the  velocity 

potential  for  a  steadily  moving  pressure  distribution  of  oscillat¬ 
ing  strength.  Lunde  [ 1951b]  has  investigated  the  special  case  when 
K#(  $X)m  fJ (fft ^"crPTy^nd  has  shown  that  as  t  •*  the  expressions 
(22.48)  and  (22.49),  after  a  change  to  moving  coordinates,  approach 
asymptotically  to  the  expressing  (21.26)  or  (21.31)  proper iy  modi¬ 
fied  for  circular  symmetry  (the  assumed  symmetry  is  not  essential). 
The  computation  is  Interesting  but  will  not  be  carried  out  here. 
This  procedure  for  obtaining  (21.26)  or  (21.31)  yields  the  velocity 
potential  without  necessitating  the  extra  boundary  condition  re¬ 
quiring  the  motion  to  vanish  as  X  +■ 00  • 

As  was  mentioned  in  connection  with  the  solution  of  the 
Cauchy-Poisson  problem,  the  Green's  function  for  some  other  simple 
configurations  can  be  found  by  the  method  of  reflection. 

The  complex  velocity  potentials  for  two-dimensional  motion 
which  correspond  to  (22.48)  and  (22.49)  are  a®  follows: 


(22.50) 
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Certain,  special  cases  have  been  investigated  in  more  detail. 
Stoker  [1953]  has  treated  the  motion  resulting  when  a  pressure  dis¬ 
tribution  constant  in  time  for  t  >0  is  suddenly  applied  to  a  uni¬ 
formly  moving  stream  of  depth  h.  The  velocity  potential  may  be 
obtained  from  (22.51)  by  taking  p(%fT)~  P*($-cT)  and  transferring 
to  moving  coordinates.  His  aim,  as  was  that  of  Lunde  in  the  com¬ 
putations  described  above,  was  to  show  that  the  potential  (21.40) 
can  be  derived  without  a  special  assumption  about  its  behavior  as 
X,  -»  +-£*o  .  The  same  can  be  carried  through  with  (22.50)  to  derive 
(21.38).  If  one  assumes  p( £,T)c  fo,(^)c4»ffTtpr(^j/^'a^then  one  may  also 
derive  (21.21)  or  (21.23)  from  (22.49)  or  (22.50),  respectively, 
as  asymptotic  expressions  for  large  t  without  having  to  impose  a 
radiation  condition.  Voit  [1.957b]  has  investigated  the  surface 
profile  for  large  t  when  f>($,T)*p(V  for  £<r  cr  <cT,  f(f,r).o 
for  ji  >  C  r-  or  for  T  >  T  . 

Waves  resulting  when  a  body  is  set  into  motion.  Many  of  the 
problems  solved  in  sections  17-20  by  means  of  source  distributions 
can  be  formulated  as  initial-value  problems  and  solved  by  the  same 


procedure  If  one  uses  the  appropriate  time -dependent  Green's  func¬ 
tion.  We  shall  consider  briefly  several  examples,  omitting  details. 

In  (19.28)  the  velocity  potential  was  given  for  the  motion 
resulting  from  an  oscillator  in  a  wall,  described  by  (19.26).  It 
•as  assumed  there  that  a  steady  situation  had  been  reached  in  which 
the  motion  was  purely  harmonic  in  the  time.  Suppose  instead  that 
the  motion  of  the  oscillator  described  by  (19.26)  is  to  start  at 
t  -  0  and  that  for  t  <f0  the  oscillator  and  fluid  are  at  rest.  It 
Is  easy  to  verify  that  the  time-dependent  velocity  H1*' 7'j'^otent ia  1 
is  still  given  by  (19.28)  if  for  the  Green's  function  G  one  uses 
(13.50)  with  m  -  1.  The  last  term  in  (13.50)  will  give  the  trans¬ 
ient;  aspects  of  the  motion.  For  two-dimensional  motion  the  time- 
dependent  wave -maker  has  been  considered  by  Kennard  [1949],  who 
also  gives  an  estimate  of  time  necessary  for  the  transient  terms 
to  die  out. 

In  (20.65)  the  velocity  potential  has  been  given  for  a  "thin" 
ship  moving  with  constant  velocity  c;  it  is  assumed  there  that  a 
steady  state  has  been  reached.  Let  us  now  suppose  the  same  ship 
to  move  with  velocity  c(t),  t  >0,  but  that  it  and  the  fluid  have 
been  at  rest  for  t  <0.  As  in  (20.64)  we  take  a  coordinate  system 
moving  with  the  ship.  Then  from  (20.26)  it  follows  that  the  velo¬ 
city  potential  ht)  must  satisfy  the  boundary  condition 

±  a,  t)  *  p  e,(t)  Fjc,  (X,  y) . 

A  Green’s  function  enabling  us  to  construct  <j>  can  be  easily  ob¬ 
tained  from  either  (13,49)  or  (13.53).  However,  let  us  take 
c(t)  -  C,  a  constant,  for  t y  0,  l.e.  we  suppose  the  ship  to  attain 
Instantaneously  its  final  velocity.  The  Green's  function  for  this 


situation  has  already  been  written  out  explicitly  in  (13.51).  Set¬ 
ting  there  <2»*  ,  A#  »  y  ,(**3  an<*  call,  ling  the  resulting  func- 
tion  the  velocity,  potential  for  the  problems  at  hand 


Is 
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(22.  52) 


Having  found  otic  may  then  compute  the  force  upon  the  ship 
and  obtain  formulae  analogous  to  <20.67)  or  (20.69).  The  computa¬ 
tions  for  infinite  depth  was  originally  made  by  Sretenskii  1 1937); 


Lunde  [ 1951a]  gives  an  exposition  of  this  result  and  extends  it  to 
include  thin  ships  moving  in  an  infinite  expanse  of  fluid  of  depth 
h  and  down  the  center  of  a  canal  of  width  b  and  depth  h.  In  these 
computations  c  is  allowed  to  be  an  arbitrary  function  of  t.  We 
refer  to  Lunde 's  paper  for  the  results. 

Havelock  [1948,  1949]  has  considered  the  accelerated  motion  of 
a  submerged  horizontal  circular  cylinder  in  fluid  of  infinite  depth. 
The  complex  velocity  potential  is  expanded  In  a  Laurent  series  about; 
the  center,  starting  with  a  dipole.  In  order  to  satisfy  the  other 
boundary  conditions,  one  makes  use  of  (13.54)  to  obtain  singulari¬ 
ties  of  the  proper  sort.  The  boundary  condition  on  the  circle  then 
yields  an  infinite  set  of  equations  for  determining  the  coefficients 
in  the  Laurent  series.  After  finding  as  many  terms  as  seems  neces¬ 
sary  for  a  suitable  approximation,  one  may  compute  the  force  on  the 
cylinder.  Havelock  has  carried  this  out  for  the  first  two  singulari¬ 
ties  (a  slight  inconsistency  in  the  approximation  is  corrected  in 
Maruo  [1957])  and  has  made  numerical  computations  for  an  Impulsive 


start  and  for  a  constant  acceleration.  Consider  an  impulsive  start 
with  instantaneous  Acceleration  to  constant  speed  c,  and  let  the 
cylinder  have  radius  a  and  its  center  Ik  submerged  to  depth  h. 

Then  the  two  leading  terms  in  the  resistance  are  given  by  R0  , 
the  steady-state  resistance  given  in  equation  (20.52),  and  by  the 
transient  term 


R. 


^  /  .  ,  jr  > 


=  1 


(22.53) 


Figure  28,  taken  from  Maruo  [1957],  shows  R>tjR.0  plotted  against 
C  t/h  for  c/ ^4  -  /. 
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FIGURE  28 

An  exposition  of  the  theory  of  accelerated  motion  of  submerged 
bodies  is  given  by  Maruo  [1957,  ch.  3].  Both  two-  and  three-dimen¬ 
sional  problems  in  fluid  of  infinite  or  finite  d<  pth  are  considered. 
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We  note  that  the  use  of  Koch in’s  H- function  may  be  extended  with  no 
difficulty  to  tine-dependent  motion;  this  has  been  done  by  Karoo 
and  earlier  by  Has kind  [1946b]. 

An  investigation  of  Palm  [1953]  also  fits  'into  the  category 
of  problems  under  consideration.  In  considering  flow  over  an  un¬ 
even  bottom  in  section  20  <x  ,  it  was  necessary  to  impose  an  upstream 
boundary  condition  in  order  to  obtain  uniqueness  if  the  velocity  is 
subcrt tical .  In  order  to  avoid  this  extra  condition  he  formulated 
an  initial-value  problem  in  which  the  fluid  is  at  rest  and  the 
bottom  suddenly  starts  to  move  The  asymptotic  expression  for 
large  t  in  a  coordinate  system  moving  with  the  bottom  agrees  with 
the  results  in  20  cy  . 

Tsunamis  r  and  submarine  explosions.  A  tsunami  Is  an  ocean  wave 
originating  from  a  sudden  upheaval  or  recession  of  the  ocean  floor. 
If  one  assumes  an  ocean  of  uniform  depth  h  and  if  the  disturbance 
occurs  in  a  region  S  of  the  bottom,  one  may  approximate  this  situa¬ 
tion  by  the  boundary-value  problem  in  which 
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If  the  time-interval  of  the  disturbance  is  short  (i.e.»  if  A- 
is  small),  the  solution  for  <f>  is  given  approximately  by  distribut¬ 
ing  over  S  sources  of  a  form  easily  derived  from  (13.53) .  In  fact, 
in  (13.53)  let  a  *  \  ,  4*'  L  ,  C  -  £  ,  and  let 
denote  the  resulting  function  by  ).  Then 
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is  the  approximate  solution.  If  one  assumes  Vfa'yt)*  W**  j )  for 
t  C  r  ,  then  </«,  takes  the  following  simple  form  for  t?T: 


14 


fj  t~"b  i  j^/ 


w'h’sr©  0"  -  <•?&  towJ*  .  If  the  deformation  la  assumed  to  take 
place  so  quickly  that  one  may  let  T-»  0  while  keeping  VT»L($,$ con 
at ant  (i.e, ,  keeping  the  same  total  deformation),  (22.56)  become* 


d^/ 

(22.57) 


and  the  solution  (22.55)  is  no  longer  approximate  for  the  formulated 
problem. 

A  further  approximation  may  be  obtained  by  assuming  the  area 
of  disturbance  to  be  so  localized  that  on©  may  assume  the  whole 
disturbance  to  originate  at  one  point,  say  Then  (22.55) 

becomes  simply  (22.57)  with  L  replaced  by  Qm$id$df\ and  £ 

Although  this  may  be  a  reasonable  approximation  to  the  explosion 
of  a  mine  on  the  ocean  floor,  it  is  not  in  general  suitable  for  a 
tsunami  since  the  diameter  of  the  region  of  disturbance  in  the 
latter  may  be  many  times  the  depth  of  fluid. 

A  comparison  of  (22. 55),  with  (22.57)  for  ^3jwith  (22.15)  shows 
that  one  may  expect  the  same  qualitative  behavior  for  tsunamis  as 
for  waves  resulting  from  an  initial  deformation  of  the  free  surface. 
In  fact,  if  one  makes  the  substitution  (22.19)  in  the  expression® 
for  the  surface  profiles,  they  become  the  following,  respectively, 
for  the  initially  displaced  surface  and  the  tsunami: 
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On©  may  study  the  development  of  7  along  the  lines  worked  out  in 
section  15  for  two  dimensions. 

Many  of  the  investigations  of  tsunamis  have  been  devoted  to 
an  examination,  of  the  profile  for  a  given  type  of  initial  bottom 
disturbance.  The  classical  papers  on  tsunamis  are  by  Sano  and 
Hmsegawa  [1915]  and  Syono[l936).  They  have  recently  been  investi¬ 
gated  by  Takahashi  [1942,  1945,  1947],  Ichiye  [1950],  Gazaryan [ 1955] 
and  others.  Since  the  shape  of  the  bottom  and  the  configuration  of 
the  shore  are  of  obvious  importance  in  a  geophysical  application  of 
the  theory,  much  recent  attention  has  been  given  to  this  aspect  of 
the  propagation  of  tsunamis.  Gn.go.rash  [  1957a]  hat*  gi,  ,/en  a  brief 
survey  of  the  literature  together  with  a  substantial,  bibliography. 

The  waves  resulting  from  an  exploding  submerged  mine  may  be 
represented  approximately  by  using  the  velocity  potential  for  a. 
source  whose  strength  «M,  t)  has  the  form  of  a  square  pulse  of 
duration  T.  One  may  then  determine  ^  from  either  ( 1*3 . 49)  or  (13.53). 
If  one  assumes.  T  very  small  and  forms  the  limit  as  TWO  while  keep¬ 
ing  mT-  Q  constant,  one  finds  easily  the  following  expressions, 
for 

infinite  depth. 

<p  *  Z(?/e4’  ^ \J.(<tR)<r  fa)  die  /  <t\  A  ; 


(22,59) 


depth  h. 
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Again  one  may  examine  the  development  of  the  surface  profile  by  the 
method*  developed  in  section  15. 

Investigations  of  waves  generated  by  a  sudden  pulse  of  the 
above  or  similar  sort  have  been  made  by  Lamb  [1913,  1922]  and 
Terazawa  [1915];  both  took  the  fluid  to  be  infinitely  deep. 
Sretenskli  [1950,  1949]  has  made  a  similar  study  when  the  source 
(two-dimensional)  is  situated  on  the  bottom  of  a  rectangular  basin 
and  within  a  fluid  layer  covering  a  solid  sphere.  Sezava  [ 1929a, b] 
has  included  the  effect  of  compressibility  of  the  fluid. 

One  should  recognize  that  such  studies  can  elucidate  only  n 
small  part  of  the  phenomena  associated  with  underwater  explosions. 
An  investigation  of  the  migration  and  oscillation  of  the  explosion 
bubble  requires  different  analytical  methods.  Furthermore,  if  the 
explosion  is  too  violent  the  linearized  boundary  condition  on  the 
free  surface  may  not  be  a  useful  approximation. 

Freely  floating  bodies.  The  motion  of  a  freely  floating  body 
following  an  Initial  displacement  is  of  considerable  interest  and 
practical  importance,  but  also  leads  to  a  difficult  mathematical 

P  # 

problem.  Uniqueness  of  solution  follows  from  the  argument  in  sec¬ 
tion  22  ex'  .  For  the  sake  of  perspicuity  let  us  restrict  ourselves 
to  motion  constrained  to  be  vertical,  i.e. ,  heaving  motion.  Then 
from  (19.59)  and  (19.62)  the  boundary  conditions  to  be  satisfied 
on  the  surface  of  the  body  in  its  equilibrium  position,  , 
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*  «  S.  (22.62) 

(The  notation  is  explained  in  section  19  ft.)  In  addition  <p  roust 
satisfy  the  free-surface  condition 


(x>  °t  jti  )  +  J  (X,  O,  J,  i)  -  o 

and  initial  conditions,  say 
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(22.64) 


(22.65) 


As  in  many  previous  cases  one  maj  reduce  the  problem  to  the 
solution  of  an  integral  equation  by  use  of  a  Green's  function.  In 
(13.49)  replace  (a,b,c)  by  ^  and  m(t)  by &($,*! t$,i)  \  denote  the 
resulting  function  by  : 
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We  now  attempt  to  express 
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by  the  integral 
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for  than  <22.63)  and  (22.64)  will  be  satisfied,  One  should  note 


especially  that  the  relation  of  j>  to  if  is  more  complicated  here 
than  in  problems  typified  by  (16.12),  for  the  past  history  of 
!•  involved  In  .  The  conditions  (22.61)  and  (22.62)  now  be¬ 
come 
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where  Jf  also  enters  into  the  equations  through  .  The  two  equa¬ 
tions  form  a  pair  of  coupled  integro-dif ferentia .1  equations  for  K 
and  u ,  .  I  t  is  evident  that  one  can  probably  not  hope  for  an 
analytic  solution  even  for  simple  configurations. 

Sretenskil  [1937b],  for  two  dimensions,  and  later  Haskind 
[1946b]  for  three  dimensions  simplified  the  problem  further  by 
assuming  the  body  to  be  ’’thin",  i.e,,  if  the  surface  is  given  by 
y~  1'  F fa ^ )  *  bj  replacing  the  bo’-ndary  condition  (22.61)  by 

±  0,{)  -  T  yi(t)Fy(\X,%)  (22.70) 

and  S.  by  the  projection  of  on  the  plane  *  -  0  [cf.  (20.26) 

and  (20.64)].  With  this  further  assumption  one  can  immediately 
satisfy  (22.68)  by  takinp 
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Equation  (22.69)  then  becomes  an  integro-dif ferentia!  equation  for 
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The  procedure  ie  open  to  eone  objection  in  that  the  substituted 
boundary  condition  (2:2.70)  does  not  sees  to  fit  into  the  femoral 
perturbation  scheme  as  developed  in  sections  10  «  ,  19  *<  and  20  ■ 

It  is  thus  not  clear  what  physical  problem  really  corresponds  to 
the  mathematical  problem.  However,  this  seems  to  be  the  closest 
■anyone  has  come  to  reducing  the  equations  to  a  manageable  form, 
Sretenskli  solved  his  resulting  integro-dif ferential  equation  nuneri 
cally  for  a  surface  described  by  ,  where 
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The  resulting  graph  of  ie  shown  in  Figure  29  with  a 

dimensionless  abscissa  £  .  In  spite  of  the  quest lonableness 

of  the  formulation  of  the  problem,  the  graph  la  instructive  in 
showing  the  difference  between  a  damped  harmonic  oscillation  and 
the  solution  of  Sretenskli ’s  integro-dif ferential  equation.  Appro*! 
mate  methods  of  solution  to  the  problem  which  assume  that  the  fluid 
motion  at  any  instant  Is  independent  of  its  past  history  lead  to 
damped  harmonic  oscillations. 
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23.  Wave#  in  basins  of  bounded  extent 

The  study  of  wave  notion  in  a  basin  presents  no  special  diffi¬ 
culties  not  already  encountered  earlier,  and  has  a  particular  In¬ 
terest  because  of  the  many  opportunities  of  observing  such  eaves. 
Certain  general  aspects  of  the  problem  nay  be  considered  as  being 
contained  in  earlier  sections.  For  example,  the  general  discussion 
of  initial-value  problems  in  section  22 cx  applies  to  motion  in  a 
basin.  However,  in  order  to  make  use  of  the  results,  in  particular 
of  equation  (22.8),  in  constructing  a  solution  one  must  have  prior 
knowledge  of  the  time -dependent  Green's  function  for  the  geometric 
boundary.  Although  the  method  of  images  can  be  used  together  with 
(13.49)  or  (13.53)  to  construct  the  Green's  function  for  certain 
simple  configurations,  an  explicit  analytic  solution  is  generally 
not  available. 

The  time -dependent  problem  has  also  been  approached  in  another 
manner  by  Hadamard  [1910,  1916],  who  derived  an  integro-dif feren- 
tial  equation  for  the  function  describing  the  free  surface. 

Hadamard 's  short  notes  have  been  worked  out  by  Boullgand  [1912, 

1926,  1927]  and  developed  further.  Certain  of  Boullgand *s  investi¬ 
gations  indicate  that  singularities  which  may  occur  at  the  inter¬ 
section  of  the  plane  0  with  the  basin  walls  are  a  result  of 
linearizing  the  free-surface  boundary  condition.  For  an  exact 
statement  one  should  consult  the  original  papers.  There  is  a  nriaf 
treatment  of  Hadamard  s  equation  in  Vergrje  [1928,  §  10,  14]. 

Moiseev  [1953]  has  developed  a  treatment  of  the  time -dependent 
problem  which  generalizes  somewhat  the  method  used  in  section  23  or . 


In  section  23  <x  we  give  some  general  theorems  concerning  Motion*# 
periodic  in  time,  and  another  solution  of  the  initial  value  problem. 

In  section  23 j3  wave  motions  for  several  special  configurations  of 
the  boundary  are  given.  In  section  23^  the  theory  of  wav©  motion 
in  movable  basins  la  considered. 

23  ck  ,  Periodic . waves  in  basins —-gene ra  1  ^  theorems . 

If  the  motion  :1s  periodic  in  time,  the  velocity  potential  may 
be  found  by  solving  a  Fredholm  integral  equation,  again  after  in¬ 
introduction  of  an  appropriate  Green’s  function.  Assume  #0* 

(f  (X  g)  c«j  (a  t  f  Tj  ;  then  f  must  satisfy  the  boundary  conditions 

f  ( X ,  0,  J  )  -  y  f(x,  O,  7  )  -  0  f  (H,  t)  F  ,  -P  -  C y2/a  , 

1  '  '  ‘  1  i;  23 .  i ) 

if 

VH,  =  0  f  (X,  j)  Ov.  S  f 

•here  F  is  the  part  of  the  plane  m «  O  occupied  by  the  free  sur¬ 
face  at  rest  and  S  Is  the  surface  of  the  basin.  Let  C  faity 
be  the  Green’s  function  for  Neumann's  problem,  i.e,, 

+  Qo} 

where  b,  is  regular  in  the  region  occupied  by  fluid  and  G  satis¬ 
fies  the  conditions 

„„  S  ,  C^x, 0,^,1, 7, 5)“  c  F  (23  2) 

where  C  is  an  arbitrary  nonzero  constant.  In  addition,  in  order  to 
make  the  definition  of  f  unique  we  require 

J ff  f  dc r  a  1 0  ■ 

S  +  F 


y 


It  then  follows  fro*»  Gross's  theorem  that 

fa')- )) -r,  5 */7 //c  rao. j;-» w 


z ; 


(23.3) 


If  one  sow  lets 


<9 


one  obtains 


ft** })-27  ]. V  o,  5) fY*  » ,  fM */, 


(23.4) 


a  homogeneous  Fredhola  integral  equation  for  rf(x,0,^)  .  If 

can  toe  found,  then  is  de te rained  by  (23.3).  From  the  theory 

of  such  integral  equations  there  will  exist  a  sequence  Vt 
eigenvalues  for  which  (23.4)  will  yield  solutions  f> ,  'fj f  ■  *.  The 
functions  corresponding  to  different  s  are  orthogonal  on 

F,  as  shown  in  (16.10).  If  several  ^  3  have  the  same  value,  the 
corresponding  *•/;  b  can  toe  orthogonal ized.  The  ft  also  form  a 
complete  set  on  F.  Each  solution  ‘f?i  yields  a  standing  wav©  in 
the  basin. 

It  is  possible  to  use  these  solutions  to  solve  the  Initial- 
value  problem  formulated  in  (22.2),  but  with  fc  m  @  ,  Let  \) 

and  7y0)  be  given.  We  trv  to  express  f(X,0,^,t)  in  the 

following  form: 


I  i-/ 


Sa*-  Cj  t 


(23.5) 


Then 

'°)~  “-2  &L  ft 

~<j  -  i>H(x  ,0,1,0) -Z  «V  O’f.'fiiw*  3) 


(23.6) 


» 


Slue®  the  for*.  it 

the  coefficient*  CLU 
Banner. 

In  order  to  use 


con  pie t©  set  of  orthogonal  functions  over  F, 
and  can  be  determined  in  the  usual 

Is  then  determined  by  (23.5)  and  (23.3). 
the  integral  equation  (23.4)  on©  must  first 


find  G»  the  Green’s  function  to  a  Neumann  problem  for  a  region 
having  a  corner  along  the  curve  of  intersection  of  the  plan© 
unci  the  basin  walls.  The  difficulty  with  the  corner  can  be  over¬ 
come  in  certain  cases.  If  the  basin  wall  intersects  the  plane 
perpendicularly,  then  the  basin  plus  its  reflection  in  the  plane 
ym  Q  has  a  boundary  without  this  corner.  If  ft *} »  is 
a  Green's  function  for  the  Neumann  problem  for  the  extended  region, 


then 


t  7- ') =  +  (23.7) 

is  a  Green's  function  for  the  original  region.  For  some  other 
special  regions  one  may  construct  a  Green's  function  by  the  method 
of  images,  even  though  the  intersection  with  the  plane  y  »  O  is 
not  perpendicular. 

As  mentioned  above,  each  f,  represents  a  standing  wave  of 
frequency  .  It  may  happen,  as  we  shall  see  presently,  that 

two  or  more  0"t  ,£  are  equal.  Let  cr  be  such  an  eigenvalue  and 
Lf  ''  and  '/  ^  iwo  of  the  corresponding  potential  functions.  By 
forming  the  standing- wav©  solution 

jc«3  O't,  A  ,  f  2  2  m  1/  (23.8) 

one  may  vary  continuously  the  position  of  the  nodal  curves,  say. 

If  n  independent  correspond  to  cr  ,  then  the  possible  nodal 
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curve*  for*  m  (at-  U~  pa  rare  ter  family  of  curve*  in  F,  With  the 


i  (0 


two  solution*  (f”'  and  one  may  aino  form  the  solution 

f  Ux  y  i  J  1 1)  -  V^'V,  7*  +  f(l)(rXt^p  w  cr  f  . 


(23.0) 


The  nodal  curve*  will  now  migrate  from  those  of  to  those  of 


(f 


(l) 


and  then  on  again  to  those  of  'f ! '  '  .  If  9s 


tfSV'/y 


and 


«/Yt> 


have  a  common  aero  at,  say,  (Jt»,  <j  0)  ,  then  a  nodal,  curve  for  <p  will 
always  pass  through  (X®,  .  Near  (x®t^)  the  waves  will  appear 

to  progress  about  !-*<>. J«)  like  spokes  moving  about  a  wheel.  There 
may,  of  course,  be  several  such  centers. 


23^3^  .  ^  Some  age 


cia.1  boundaries. 


It  is  possible  to  give  explicit  solutions  for  standing  waves 
for  several  particular  configurations  of  the  basin.  The  variety 
of  such  configurations,  however,  is  rather  small.  As  a  preliminary 
we  note  that  if  the  basin  has  a  flat  bottom  at  depth  h  and  If  the 


side  walls  form  a  vertical  cylinder  making  a  section  C  with  ^  » 
then,  from  section  13  <xf  ,  we  have 

)m  (f(0C^)o»yC  C*o(crt  +  T)  f 

In,  'Hv.  y  •  O  , 

where  fVjt.j)  s  a  solution  of 

"t  X  X  +  +  lH  •  /  *  0 

satisfying 

c 


O 


(23. 10) 


(23. 11) 


t,- 0 


<23. 12) 


W 


1 


iTi 

The  boundary  condition  (23.12)  will  limit  ,  and  hence  <j  ,  to 
a  discrete  sequence  of  eigenvalues 

♦Hr*  ,  AW  J  \  t  <y ,  ,  <J\  ,  •  ■  •  (23.  13) 

In  a  coordinate  system  in  which  (23.11)  can  be  separated  it  is 
usually  possible  to  find  the  standing  waves  in  basins  whose  side 
walls  are  constant-coordinate  surfaces.  These  statements  will  be 
illustrated  below  for  rectangular  and  cylindrical  coordinates. 

In  connection  with  the  special  cases  treated  below  we  call 
attention  to  papers  by  Honda  and  Matsushita  [1913]  and  Sasaki  [1914]. 
The  authors  investigated  experimentally  in  a  systematic  way  the 
various  modes  of  motion  in  rectangular,  triangular,  circular,  ring- 
shaped,  circular-sectorii i  and  ring-sectorial  basins  and  compared 
measured  with  calculated  periods.  In  most  cases  the  agreement  is 
with  2%.  Photographs  showing  the  various  modes  were  obtained  by 
sprinkling  the  surface  with  aluminum  powder  and  exposing  a  photo¬ 
graphic  plate  for  about  one  period.  The  nodes  then  show  up  as 
dots,  the  rest  as  streaks.  In  connection  with  a  study  of  the  ex¬ 
citations  of  waves  in  a  port,  McKnown  [1953]  has  also  investigated 
experimentally  the  standing  waves  in  circular  and  square  basins; 
some  striking  photographs  are  included.  Apte  [1957]  has  studied 
further  the  theory  of  the  excitation  of  standing  waves  in  a  square 
basin  and  has  also  given  experimental  results.  Perhaps  the  first 
theoretical  investigation  was  by  Rayleigh  [1876,  pp.  272-9];  he 
compared  his  predicted  frequencies  with  observations  of  his  own 
and  of  Guthrie  [1.875] 

Rectangular  basin.  Let  the  basin  walls  be  given  by 

o,  *-*■  ,  v  1 


Them  fro#.  (13.6)  one  may  write  down  immediately  the  solution 


<P  *  A  a *L  «K»(  y  +•  t) 


C#a 


T 


Jt 


fo  V 

T 


cm(<rt  +r) 


(23. 14) 


^  ( ~J\  +■  ^  ^  *  4*v»  ♦Ho  4^,  p,  £  •  <?f  •  . 

Tims  the  choice  of  the  integers  p  and  £  determines  4h»  and  then  cr  . 
If  the  basin  is  square,  i .  e.  ,  /«  4  ,  then  the  same  values  of  **» 
and  or  may  correspond  to  two  different  solutions  obtained  by  inter¬ 
changing  f»  and  £  ,  assuming  p  f*£  .  However,  this  may  also  occur 
for  other  rectangular  basins  If  4  and  l  are  commensurate. 

Circular-cylinder  basin.  Let  the  basin  have  radius  a.  Then 
from  (13. 8)  we  find  the  solutions 


<f>  ~  A  6#»'W  4K. if  *j  4  „  R.)  C*i  (4v  «f  4  f)c+a(<rt  *T)  ,  4%  *  o,  I, 


'* . * 


t 

'Xu,  (4U.  O')  -  O  f  £  -  4*.  «  UM'M  4U  c  L  .  (  23 . 1 5) 


Thus  /Hv0  must  be  selected  so  that  4U*  <x  is  one  of  the  zeros  of  ^  ; 

this  then  determines  O'  .  Tor  n»o  the  wave  crests  and  nodes  Lie 
on  concentric  circles,  the  number  of  such  nodal  circles  depending 
upon  which  zero  of  tJ 0  is  used  to  determine  .  If  -n  >  1  , 

then  to  the  same  cr  there  correspond  two  independent  solutions 
(  <f*  o  ,  S  *  |  JT  ),  and  the  remarks  made  in  connection  with  (23.8) 
and  (23.9)  apply. 

The  standing  waves  in  a  basin  shaped  like  a  sector  of  a  circle 
may  be  obtained  from  (13.8).  If  tX#  is  the  angle  of  the  sector, 


then 


®  A  Ce>4.  4U»  ( u  f  4)  J^jr  (* %.)  Coo  cea( e rtf  l)/  <K  * 

*  "oT#  ° 

I  v 

tJ^IT  (ma,o  Ol)  -  O  /  -?  *  4H.  4  . 

a*..  ? 


>  (/'  1  V 


(23  15) 


3  bo 


If  the  basin  is  ring-shaped,  with  inn*r  radius  b  and  outer 
radius  a,  then  from  (13,8)  one  find®  (ef.  Sara®  [1813],  Campbell 
[ 1933]) : 

f*A  c**l  **o(*]+l)!  xl *2 h  J* ♦  V  Y+  (**.*)]<*»  fa*  +s)e*>(<rt+*)f 

(33, 16) 

AA**  «  d  [cu-4  •»  A,  f  4*  k  ^ji  Jf  "  ‘  "  f 

Formulas  for  sectors  of  a  ring  may  be  obtained  and  are  similar  to 

(iO.lfl). 

Basins  with  sloping  side-walls.  There  are  very  fee  explicit 
solution®  known  when  the  sides  are  not  vertical.  If  the  basin  Is 
a  horizontal  cylinder  bounded  at  either  end  by  vertical  walls  at, 
say,  \*o  and  ,  the  theory  of  progressive  waves  in  canals, 

developed  in  section  18  jf  ,  can  be  carried  over  with  only  small 
changes,  namely  replacement  of  4j  cvafWfr) where  now 

it  is  restricted  to  the  values  and  cr  correspondingly.  Thus 

( 18.39)  and  (18.43)  give  the  velocity  potentials,  after  the  indi¬ 
cated  modifications,  for  various  modes  of  oscillation  of  a  fluid 
In  a  basin  of  triangular  section  whose  sides  form,  an  angle  of  45° 
with  the  horizontal.  However,  even  though  these  formulas  may  be 
used  also  for  the  two-dimensional  modes,  when  4*0  ,  they  do  not 
give  the  gravest  two-dimensional  mode  except  by  a  limiting  process 
(described,  eg.,  in  Lamb  [1932,  p.  443]). 

The  two-dimensional  modes  of  motion  in  triangular  basins  whose 
sides  form  an  angle  with  the  horizontal  may  also  be  studied 

by  use  of  the  methods  introduced  in  section  IT  £>  for  standing  waves 
on  beaches.  Indeed,  It  is  apparent  that  Kirchhof f  [1879]  considered 
his  investigation  of  waves  on  beaches  as  a  preliminary  to  the 
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problem  at  hand.  Because  his  approach  la  systematic  we  shall  des¬ 
cribe  it. 


In  order  to  use  the  results  of  17 ^  we  take  one  side  as 


<r 


•  V  I 


I 

t 


3 


4  C  *  ^  ;  let 


the  other  side  be  given  by 


2a  **  Jf  t 


*>* 


(23. 17) 


Then  the  complex  potential  -ftp  must  satisfy  not  only  (17.31)  and 
(17.32),  but  also 

-f(2as  ~  A  t  ^ )  *  {  (trCk  -  JlC  '  )  t  (23. 18) 

which,  taken  with  (17.31),  implies  that 


-ffjp-  f  (j 


’i‘1t  *Za.rlt( l-eU*)) 


(23. 19) 


In  order  to  satisfy  (17.31),  (17.32)  and  (23.19)  Kirchhoff 
first  takes 

(23ao) 

Substitution  in  (17.32)  yields  (with  fl  *exf>-2tf  as  before) 


!-  fl'  =  O  .  3. 


.  n  tn. 

tV  I  -t/S  ,  ,  y. 

. t  Ju  *  i  3.  .  /  +  /- 


4/+  *4 


I  -~W  1  | 


(23.21) 


Thus,  since  ^-4^1"/^,  on©  must  have  4/  "  'Hy  p  * ; and  -  a 

an  integer.  If  one  takes  * o  ,  then  (23.20)  becomes 


4 


<-xi 


4  (jf 


(23.22) 


Condition  (17.31)  requires  B0  to  be  real.  Condition  (23.18)  or 
U3.19)  remains  to  be  satisfied.  The  function  fo)  In  Its  assumed 


form,  is  apparently  overdeternined,  and  it  is  possible  to  show 
that  for  j?  ?Z  not  all  conditions  can  be  satisfied.  For  g M  2,  *Kml 
and  <?•  3^,  Av.  *  1  ,  (23.18)  can  be  satisfied.  The  potential  functions 
are  as  follows: 


X  •  fi/4 

■  £.[(/-  vx)0  +  vy)  -  iv  (x  +  y)  ( t-  1  »fx-%))] / 


(23.23) 


!/v/  *  *h 


/ 


:  (23.24) 

>- 1 ,  t.  [2  * i  -  r?  +1)1] 

.  (dx  .  v*fi)(vx  -  VuTi  ■■  lfi)(vx-{l)] ,  O -  m  /"i f  4  *  , 

Here  h  is  the  depth  of  fluid  at  the  deepest  point.  The  surface 
profile  for  ^  -  45°  is  a  straight  line,  for  ^  -  30°  a  parabola. 

In  order  to  find  the  higher  modes  of  oscillation  Kirchhoff 
returns  to  the  form  (17.33)  for  ■  It  then  follows  as  before 

that  (17.34)  must  hold  and  that  r  must  be  even,  say  2g  .  Now, 
however,  instead  of  taking  A  =  1  it  is  left  to  be  determined  by 
(23.19).  Substitution  of  (17.33)  into  (23.19)  gives 

L*. + 1 

^(t,+2  ~  ^4  ^  A  va(!>  ( 1- a)  1  h.  -  Of  1 ,  •  *  *  ,  M  -  5  (23.  25) 

Altogether  there  are  then  - 2 -2-n - 5  independent  equations  to 

determine  A,  ,  .  and  also  A  and  va,  Again  the  conditions  can 

be  satisfied  for^-JT/4  and  J/c  ■ 


Th«  solutions  are  a«  follow,  shore  C  Is  as  arbi 

trary  real  constant 


{ /§ )  -  C [ Cea  A  v  ( j  ♦  6*->A  A  v  ()  V)  ] 

A  *  cWi,  A  vft  »  -  o*4  An)G*  /  (33.33, 

f/'j )  «  (f  f"'c#o  A  '»  ()  -  &  (I  ”  fc  >)  "  o**^  ^  J 

/I  *  J«maX  A  va  *  A  v cu  , 

The  values  of  A  and  u  can  easily  be  date  mined  graphically.  For 
the  first  set  of  solutions  the  values  ofAva^wlll  be  slightly  sors 


than  5JT/4  ,  lWf4  ,  ,  .  ,  ,  for  the  second  set  slightly  less  than 


7 jr/4  ,  ft  JT/4  . .  .  These  two  sets  of  solutions  correspond, 

respectively  to  (18.39)  and  (18.43)  with  -4*0  ;  the  eigenvalues 
A  may  he  identified  with  and  n(/v  ,  respectively.  Klrchhoff 

and  Hanseuano  (1880)  carried  out  an  experimental  Investigation  of 
the  first  three  antisymmetric  modes  (equation  (23.23)  gives  the 
first  one);  they  compare  frequencies  and  positions  of  maxima  and 
minima.  The  agreement  seems  satisfactory,  although  correction® 
for  surface  tension  were  necessary  for  the  two  higher  inodes. 

The  solution  for  ^  -  30°  is  the  following: 


/fj)-  cbxr 


,  IA  I 

_  £  V  — - g  " 


(23.27) 
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.t/b  Ay[  %-<*■- fi  **■] 
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A  +  /  w  A  -  / 

T  mi  f*-! 

where  C  Is  an  arbitrary  real  constant,  /l  =1  (f~  i  (*,),  ft,  » 
and  the  eigenvalues  for  A  and  v  are  determined  by  the  equations 


All, 


■j&CWAva, 


t(}+/i>)v4fi>AvA' 


A. 


■ » -f  t,  a  >^2^ 


If  A  ia  a  solution  of  I 23.26* ,  then  also  ~  A  A  ,  y6A  and  ft  A  are 
solutions.  There  exists  a  real  solution  which  may  be  found  from 
the  equations 

,  r—  ,  „  ,  •>  i*<vj u.  iT5  A  va  **  Vj5»  jiAujA  . 

The  other  solutions  which  may  be  generated  from  these  do  not  lead 
to  expressions  different  from  '23.27'  .  The  first  eigenvalue  for 
Ava  is  a  trifle  to  *he  right  of  /Z  .  The  form  of  the  free 
surface  corresponding  to  (23.27)  is  given  by 


&*>•(  '{&*.*»(*-*■>* 


(,v 


J  (3  Avx 


A  - 1 
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(23.30) 


f*,  A  v(x  -  Za) 


,  A  +  if 


+  X-7‘ 


i  h  ^  j  ^  (A  j.ji  j  Cert  * T  ). 


Note  that  (23.24)  and  (23.27)  both  give  only  symmetric  modes. 
Macdonald  [  1896)  states  that  antisymmetric  modes,  if  they  exist, 
cannot  be  represented  in  the  assumed  forms  (23.20)  or  (17.33). 

Vint  [1923)  has  succeeded  in  finding  an  infinite  number  of 
modes  of  motion  in  an  inverted  four -sided  pyramid,  each  of  whose 
sides  makes  a  45°  angle  with  the  horizontal.  We  refer  to  the 
original  paper  for  the  exact  formulas. 

Additioral  solutions  have  beer,  obtained  by  inverse  methods 
by  Sen  [1927)  and  by  Storchi  [1919,  1952].  Storchi's  result,  al¬ 

though  restricted  to  two-dimensional  motion,  is  neat.  Suppose 
that  the  form  of  the  free  surface  is  given  as  *j(jc,  t)  -  */<'•# > -/<T  * t  r) 

_  f'(x) -n*v(crt-f  T)  ,  where  f(A')  is  analytic.  Then,  since  yfr)** ?  $(*>$ 
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f  u-co)*  %(*,<>)-'  f/j£,o>s  ^f(x,o)^[pttfx)-^F  <JO] 

and 

-<  a)  *  ™  [r‘(x )  -  ^  F/>  }]t  c*K*t 

We  may  take  the  constant  as  zero  and  write 


{(*+<])  - 

there  FM  Is  the  analytic  function  dotenolned  by  fW 
this  we  have 


(23.31) 


From 


fc,3)>z4  j)-ip[F0*ij)  +  n*-‘jXl]f 


(23.32) 


Any  streamline,  defined  by  ~t  -  real  const.,  can  now  be  taken  as 
determining  a  possible  basin  shape  corresponding  to  the  assumed 
standing-wave  profile.  Storchi  applies  the  procedure  to  several 
special  choices  of  F.  An  obvious}  disadvantage  of  this  method,  as 
well  as  of  Sen’s,  is  that  only  one  mode  of  motion  is  obtained  for 
a  resulting  basin  shape. 


23  |  .  Waves  in  movable  basins. 

In  several  preceding  sections,  especially  19  and  22  (f  ,  se 
considered  the  save  notion  occurring  In  the  presence  of  an  oscll- 
lating  body  when  the  fluid  is  exterior  to  the  body.  One  may 
attempt  analogous  problems  when  the  fluid  is  situated  inside  the 
body  Such  problems  occur  in  many  situations  of  practical  interest 
for  example,  the  sloshing  of  oil  in  ?  partly  filled  compartment  of 
a  tanker  and  the  sloshing  of  fuel  in  an  airplane  or  rocket.  In 
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each  of  these  cases  Interest  centers  upon  the  dynamics  of  the  whole 
system  as  well  as  upon  the  effect  upon  the  walls  of  the  container. 

A  further  interest  in  such  pronlems  arises  from  the  interpretation 
of  the  experiments  on  standing  waves,  referred  to  earlier,  carried 
out  by  Honda  and  Matsushita  [1913],  Sasaki  [1914],  and  Kirchhoff 
and  llansemann  [1880].  The  results  were  intended  for  comparison 
with  theoretical  prediction  of  standing  waves  in  fixed  basins. 

The  waves  were  actually  generated  by  oscillating  the  basin  and 
finding  the  frequencies  at  which  resonance  appeared  to  occur. 

We  shall  not  consider  the  most  general  motions  of  the  basin 
consistent  with  linearization  of  the  free  surface  conditions,  but 
shall  limit  ourselves  here  to  a  particular  problem  with  small  hori¬ 
zontal  oscillations.  In  section  26  of  small  vertical  oscillations 
of  the  basin  will  be  considered.  The  general  problem  of  motion 
of  a  body  containing  fluid  with  a  free  surface  has  been  treated  by 
Moiseev  [1953]  and  Narimanov  [1956,  1957],  However,  both  are  pri¬ 

marily  concerned  with  small  oscillations,  Krein  and  Moiseev  [1957] 
have  also  considered  certain  mathematical  aspects  of  this  problem. 
Okhotsimskil  [1957]  and  Rabinovich  [1957]  have  considered  the 
special  case  when  the  fluid  is  situated  In  a  vertical,  or  almost 
vertical  cylinder;  Narimanov  also  gives  special  attention  to  this 

case.  (Publication  of  the  work  of  these  three  authors  was  apparently 

is 

delayed:  it, stated  that,  for  the  most  part;,  it  was  carried  out  in¬ 
dependently  of  and  prior  to  Moiseev’s  papers.  A  particular  prob¬ 
lem,  the  ore  discussed  below,  was  treated  by  Sretenskii  [1051]  and 
later  by  Moiseev  [ 1952a, b,  1953 j.  Two  later  papers  by  Moiseev 
[ 1954a, b]  apply  the  theory  to  engineering  problems  especially 


3  (si  7 

ahips.  Waves  resulting  from  a  special  type  of  forced  oscillation 
©*  *  rectangular  tank  have  been  studied  by  Binnke  ( 1941]  nnd 
Tamiya  [1958]. 

Waves  in  a  baa  In  with  elastic  rjato^rlng  force,  Consider  the 
configuration  shown  in  Figure  30.  The  coordinate 
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FIGURE  30 

system  0  YY  is  fixed,  the  system  0 X  Y  moves  with  the  carriage. 

Let  X,  (t )  -  Go,  Ki,»  m  Xe  .  The  bottom  of  the  fluid  is  at  ^  ^  , 

the  side  walls  at  X  =  ±  Qy  .  The  motion  will  be  takers  as  two- 
dimensional.  Denote  the  mass  of  the  carr i«ge ,  per  unit  width,  by 
**€  »  that  of  the  fluid  by  4*^  and  the  total  by  rfw  «  4we  ♦  . 

Let  the  spring  conn t ant  be  /w%,  *&  ,  We  suppose  as  usual  that  the 

motion  may  be  described  by  a  velocity  potential  i).  Follow¬ 

ing  the  notation  at  the  end  of  section  2  let  $(SL,y,t)  describe 
the  motion  relative  to  the  basin,  i.e. 

|  t  )  *  $  ( 'Yf  ^  ft  )  +  14  a  X 


(23.33) 


X  lx 


<11 


1  6  ? 


le  shall  assume  that  X0  and  M0  art  tooth  small,  and  of  the 
same  order  as  <p  ,  i.e.  ,  in  the  notation  of  section  10  ec  ,  we 
assume  expansions  of  the  form 


f  i ) 


(•) 

f  Aa„  •  £  Ut  , 

f  o >  i  r(i) 

+■  £ 


(23.34) 


X*  *  £  X, 

^ "  --  It.  <j>  *  +•  £  1  + 

We  omit  the  formal  details  of  substitution  of  the  perturbation 
series  in  the  exact  boundary  conditions.  They  lead  to  the  follow 
lug  linearised  boundary  conditions  for  S  ; 

on  _a 

4),  (*■ ,  0,i)  +  cj  ,  o,t)  t  U.  x  -  0 , 

4»A  (±  tX'ity't)  s.  O, 

fg  (JC  r  K  * )  *  0  ■ 

The  pressure ,  after  discarding  higher-order  terms,  is  given  by 

(^36) 


(23.35) 


The  motion  of  the  carriage  is  determined  by  the  equation 


,  ,  y*  % 

I'W  ^  JC  q  Bt  J  jl  )  ot  —  'Vi'v  ■Cfj  X  Of  (23.37) 

where  the  integral  is  taken  over  the  wetted  surface  whi.  a  the  system 
is  at  rest.  Substitution  of  (23.36)  gives 


(23.38) 


(Equation  (23.38)  Is  also  a  direct  consequence  of  conservation  of 


momentum. ) 


AL 


i  <4*  II  4>n 


must  be 


The  velocity  potent it  1  ^  end  the  displacement  X® 
determined  together  from  equations  (23.35)  and  (23.38)  and  either 
initial  conditions  or  the  further  assumption  that  the  motion  is 
harmonic  in  t. 

As  a  preliminary  we  shall  first  suppose  that  the  basin  motion, 
i.e,  X„  ,  is  given,  so  that  only  (23.35)  need  be  satisfied.  One 
ma>  try  separation  of  variables  and  express  in  the  form 

I  '  £  %.(*)  XJJc)  (23.39) 

Laplace’s  equation  and  the  last  two  condition  of  (23.35)  are 
s at. is f lei  by 


X  Y  ~  vx-  kXr  m 

2  a,  2  2 ■  A-  ^  (23 

X,  Y  -  ** 

1  Kc  I  & 

In  order  to  find  the  corresponding  I  ,  expand  x  in  a  Fourier 


series : 


m£> 

x  Xto  (z*~+nxr 


2-  ¥  1 

l  a. 


jr  jl 


(23.41) 


and  substitute  (23.39)  and  (23.41)  in  the  first  condition  of  (23.35) 


Z  [t  1  xrX  L  J  T21 

,  j-  2m-  Z <5u  J  uOs  t&s  J  2  <a> 


(23.42) 


I  „/ 


-p±  jrL  r  T  .  O  31  .rtv  -+-x.(-»I— J 

2  to  J  tn  £  <5t-  J  i  <?<» 


0 


Let  us  set 


=■  H  r  la  i  1  n  Zcu 


(23.43) 
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Then  equation  (23.42)  yields  the  infinite  »et  of  differential  ©qua 
tions 


Ti~-0> 


( 23  *  '4  '4 ) 


T, 


2«.fi  +l|3Z^+.(  ^H  +  l  *0 


The  solution  of  the  first  set,  *'  ^v)»  i®  independent 

of  the  notion  of  the  basin  and  yields  the  symmetric  inodes  of  oscil¬ 
lation  in  a  fixed  basin.  The  solution  to  the  second  set  may  also 
be  found  by  elementary  methods,  but  will  not  be  given  here.  How¬ 
ever,  w©  note  that,  if  X0  is  harmonic,  it  confirms  that  resonance 
occurs  at  the  frequencies  of  the  asymmetric  modes  for  a  fixed  basin, 
We  now  turn  to  the  joint,  solution  of  (23.35)  and  (23.38). 
Substitute  (23.39)  into  (23.38).  Then,  after  evaluating  the  inte¬ 
gral,  one  finds 


■vw  X 


±Q  '  7  -~l)V  sw4  H  k  t  +  M*  k'x*  -  o  . 


(23.45) 


The  equations  (23.44)  and  (23.45)  taken  together  may  now  be  used 
to  determine  Jd,  and  the  .  If  we  formulate  an  initial-value 

problem  by  requiring,  say, 

_  JT 

Jt,(o)*C,,  X.(o)  =  ot  ~Vo)^o  ,  ^  (*,  ° (23.46) 

then  the  7^  are  all  zero  and  the  Tz  ^  and  must  be  devermined 

from  the  differential  equations.  As  usual,  on©  looks  for  a  solu¬ 
tion  in  the  form 


“  I  <A»  t 

Ji.  *  c  t 


'-r—  j 

1  M.  *+  I  ~  ^1^*1 


•  4  CO 


(23.47) 


•here  both  C  and  <di may,  of  course,  be  comple*.  Substitution 
In  (23.44)  and  (23.45),  followed  by  elimination  of  +  ,  , 

yields  the  following  equation  for  determining  GO  : 


4l.  il 


*  ' 

z... 


£  1  IK.  C  |  )  * 


(  <kmX\  *-1  Jr  L 


(23.48) 


One  may  find  the  solutions  graphically  by  plotting  each  side  of  the 
equation  as  function  of  w"  .  Figure  31  gives  a  qualitative  idea 


FIGURE  31 


Ill 


of  the  dim  trihut  Ion  of  solutions  *o.,  <*}, .  As  v  -*  , 

t  i, 

ejj ^ (  -  er  t f ,  .♦  0' i  this  fact,  which  can  bt  proved  analytically  and 
which  seems  clear  from  Figure  31,  would  not  have  been  so  evident 
If  we  had  divided  (23.48)  by  oJ4  before  plotting.  A  point  of  im¬ 
portance  is  that  there  is  no  intersection  for  <*->  <  °  ;  as  a  result 
the  motion  is  stable.  This  may  be  proved  as  follows.  Since  x  Vanli*  « t, 
the  right  hand  side  of  (23.48),  for  uS <  o  ,  is  greater  than  or 
equal  to 


5  2  q?  f 
77  3  *v\ 


u  /  _ •  _ . 

')l 


•O 

<r 


2.  /U  •*-/ 

Z-  a. 


F  l  ~ 


4w  x  „  x 
- -i  <0  >  LO 


(23.49) 


Hence  the  line  to1-  it  '  lies  below  the  left-hand  branch  of  the  curve 
for  it.1  ^  0  .  The  eigenvalues  to0  depend  upon  the  parameter*  4t.  , 

la.  /(^  and  2  ■**  .  Figure  32  from  Moiseev  [1953]  shows 

the  dependence  of  the  fundamental  mode  uJ*  upon  l^/ti  for  two  values 
of  / /vw  and  {t  1  x  'l 


313 


Thu  general  solution  lor  X«  and  Tt  Is 


x,(t)*  Vt  C$ 

■J  m  6 


l  t 


4*  a 


(23.30) 


Thu  solution  of  the  initial-value  proble*  formulated  In  (23.46)  sill 
not  be  completed.  It  Involves  solution  of  infinite  sets  of  linear 
equations,  approximate  solutions  can  be  obtained  by  considering 
only  a  finite  number  of  equation*  and  variables. 

The  general  theory  of  stability  of  such  systems  is  discussed 
in  Moiseev's  1053  paper.  In  an  earlier  papers  (1952b)  he  studies 
the  special  case  of  a  basin  containing  fluid  and  serving  as  the 
bob  of  a  pendulum.  If  the  suspension  is  by  a  parallelogram,  link¬ 
age,  so  that  the  container  moves  parallel  to  itself,  the  motion 
is  always  stable;  if  the  suspension  is  by  a  rod  rigidly  attached 
to  the  container,  the  motion  may  be,  under  certain  circumstances, 
unstable. 

The  last  cited  paper  by  Moiseev  describes  briefly  the  results 
of  an  experiment  with  a  pendulum;  the  measured  and  computed  funda¬ 
mental  frequencies  for  the  two  systems  of  suspension  agreed  with 
0 . 1%. 


24.  Gravity  waves  in  the  presence  of  surface  tension. 
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Apparently  the  first  one  to  investigate  the  theory  of  waves  in 
a  fluid  acted  upon  by  both  gravity  and  surface  tension  was  Kelvin 
[  1871a, to] .  However,  many  of  the  essential  features  had  been  dis¬ 
covered  earlier  through  observation  by  Russell  [1644]  and  others; 
references  nay  be  found  In  Kelvin's  papers.  A  good  account  of  the 
classical  researches  of  Kelvin  and  Rayleigh  may  be  found  in  Lamb 

[1932,  |  265-272].  Also,  Chapter  XX  of  Rayleigh's  Theory  of . Sound 

(Cambridge,  1929;  Dover,  N.Y. ,  1945]  contains  an  exposition  of  many 
of  his  own  fundamental  researches  on  surface-tension  phenomena, 
including:  waves. 

The  chief  mathematical  complication  added  by  the  action  of  sur¬ 
face  tension  is  a  somewhat  more  elaborate  dynamical  boundary  condi¬ 
tion  at  an  interface  or  free  surface.  The  difference  of  primary 
physical  interest  lies  in  the  existence  of  a  minimum  wave  velocity 
and  of  two  wave  lengths  with  the  same  velocity.  Many  of  the  special 
problems  considered  in  preceding  sections  can  also  be  solved,  when 
surface  tension  is  acting.  However,  there  has  been  little  motiva¬ 
tion  for  carrying  through  such  analyses  for  wave  motion  associated 
with  solid  boundaries,  since  it  has  been  recognized  that  the  addi¬ 
tional  forces  would  be  small.  A  further  difficulty  also  appears 
when  the  solid  boundary  pierces  the  surface,  for  an  additional  boun¬ 
dary  condition  is  required  at  the  Intersection.  As  a  result,  most 
of  the  investigations  have  dealt  with  waves  analogous  to  those 
considered  in.  sections  14  and  S ,  15,  and  22  /*  .  In  fact, 

the  complex  velocity  potentials  for  the  Cauchy -Poisson  initial-value 
problem  including  the  effect  of  surface  tension,  has  already  been 
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given  in  equations  (22.45;  am!  (22.45).  A  topic  of  particular  geo¬ 
physical  interest,  the  stability  of  an  interface ,  will  be  dealt 
Kith  in  section  26.  Waves  on  the  surface  of  a  viscous  fluid,  in¬ 
cluding  surface  tension,  are  considered  in  section  25. 

Boundary  conditions.  The  linearised  conditions  to  be  satisfied 
at  the  interface  of  two  fluids  have  already  been  given  in  equations 
(10,7)  and  <10.8)  (we  recall  that  subscript  f  refers  to  the  lower 
fluid).  If  one  eliminates  7  from  these  two  aquations  and  makes 
use  of  the  fact  that  Laplace’s  equation  is  satisfied  on  -ther  side 
of  the  boundary,  one  has  the  following  condition. 


A^(  -  O'  -for  y  <•  0  ,  AfZ  «0  \or 


(24. 1) 


7 1 '*>  5 ) “  f, j ( x ,o,$ , t )  -  <pz y(*-o.$.+), 


(24.2) 


(24.3) 

If  the  upper  fluid  is.  absent,  one  sets  f>2  and  <pz  equal,  to  zero 
and  may,  of  course,  drop  the  subscript. 

If  the  motion  is  two-dimensional  one  may  introduce  a  stream 
function  f  and  a  complex  potential  F().t)  *  +  +  '’?’  and  express  (24.2) 
and  (24.3)  as  follows: 

%(X)  *  I  »n  F  <*o)  *  U  ^  (X  +  1  °)  (24.4) 


. . .  it  mu  inn i 


If  the  upper  fluid  Is  absent  and  If  the  lower  fluid  in  Inf i 


nitely  deep,  one  may  extend  the  reasoning  which  led  up  to  Levi- 
Civita’s  differential  equation  (22.25)  to  derive  the  following  one 
which  must  be  satisfied  for  all  ?  : 


Fh  ( J  i  t )  +  l  cj  F 


7  r  .» 

rF 


(24.6) 


Furthermore,  If  the  fluid  1®  of  constant  depth  U  ,  Clsotti's  equa¬ 
tion  (22.30)  may  also  be  extended  to  include  the  effect  of  surface 


tension ; 

Ftt(}+Lkii)+Ftt(yiU)+il[F'tyl^)~F'(yiii )]-  <  flp  "fyL  L)  -  f  Vj  -&)  ]*o 


or 


-  2  L  <  y  <  o . 


(24.7) 


Elementary  solutions .  Let  us  suppose  first  that  only  one  fluid 

is  present,  and  In  addition  that 

f  (X,  p'jii)  “  ))  C*3  (<rt  *-  T ) . 

Then  (f  must  be  a  potential  function  satisfying 

-  <r -fU'O, })t  ?  fj  +f  (24.8) 


Just  as  in  section  13  o<'  ,  we  may  separate  out  the  y-variable  and 
obtain  the  following  expressions; 
infinite  depth 


where 


and 


Y( <* 'hi"7) )  -  A e  4  f '(x,  -j ) 


/i  t.  cp  f  K '  ip  *  o 


■»  T  b 

a  *  ^  4w.  f 


(24.9) 
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depth  h: 


a  r(*,p,  <24.  io) 

Zlt  (f  f  -Vu..  Cp  «  O 

CT  m  4H »  +  -p  4h  c  )  Aw  *kt  . 

One  may  also  with  no  difficulty  construct  solutions  analogous  to 
(13.3)  and  (13.4),  namely 


where 

and 


(f(^t  *rj)  ~  A  ft  "Py  A*y )  C&> 


+  f 


(24.11) 


and 

1  ('  •A  tji  y)  «  j\  t&z  urn;  (‘-  rk)  !~f  CJt/T  J 

(24. 12) 

for  infinite  and  finite  depth,  respectively,  where  tn i  in  (24.12) 
must  satisfy 

<T%  m  (~  +•  -jy  Aia,-  )  Ax-i-v  44^;  /v 

and  ffJC,  j)  must  b©  a  solution  of 

.  J>  t  j? 

^2  4vv  7  5S 

Unfortunately,  the  set  of  function 

j  O&isfn,  4V.  •  /  4^  f  k  }  j  C&ti  AH.  t  (  ^  t  k-  y  j 
is  no  longer  orthogonal  in  general,  so  that  the  convenience  of 

general  solutions  like  (16.3)  is  lost. 


It  is  not  necessary  to  repeat  the  computations  of  section  13 
since  they  remain  unaltered.  Essentially  the  only  change  is  in  the 
relation  between  the  frequency  cr  and  the  wave  number  ^ .  Here 
the  fact  of  predominant  physical  interest  is  that  for  small  values 
of  aw,  the  relation  is  controlled  chiefly  by  the  gravitational  con- 


stint  <7  and  for  large  value*  of  «*v  by  7  ff  , 

4 

If  one  forms  two-dimensional  progressive  waves  by  superposing 
the  standing  -wave  solutions  obtained  from  (24.9)  and  (24.10),  a 
further  significant  physical  fact  appears:  the  wave  velocity 
now  has  a  minimum  for  some  value  of  /w\  >  o  ,  except  for  very  shallow 
depth.  These  facts  are  displayed  graphically  in  Figure  12  and 
further  Information  is  given  in  the  associated  discussion  (the  cur¬ 
ves  were  computed  for  water  at  20°C  and  It  m  ■  <*»  or  j  cm).  Formulas 
for  the  position,  of  the  minimum  and  various  associated  values  are 
given  for  infinite  depth  in  the  following  table;  the  numerical 
values  are  for  water  at  20°C  (T  *  72.8  dynes/cm,  ^  m  .998  gm / cur ) : 

->» u  =  b.tc 

n  l  r/rv  =  Hi  t>,- 

4, - - 

*  V  -  lb  I  Caxa/'uc 

(24-13> 

When  /v^:  there  is  no  longer  a  minimum  value  of  C  for 

/^x  y  o  ;  in  this  case  C  increases  monotonically  with  /wv  .  The 
critical  depth  for  water  is  about  .33  cm.  Except  in  this  latter- 
case  every  value  of  c  has  associated  with  it  waves  of  two  different 
lengths,  each  of  which  travels  with  velocity  C.  Kelvin  suggested 
that  the  shorter  waves,  whose  behavior  is  controlled  chiefly  by 
surface  tension  be  called  "ripples".  The  suggestion  has  been 
followed  for  the  most  part  (French  "rideaux” ;  German:  "Rippeln" 
or  "Krause Iwe lien” ;  Russian;  "ryaby"],  but  they  are  frequently  also 
called  "capillary  waves’’  in  contrast  with  "gravity  waves". 
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The  relation  between  0"  and  ^  was  subjected  to  a  rather 
thorough  experimental  investigation  by  Matthlessen  [1889].  He  made 
measurements  with  water,  mercury,  alcohol,  ethyl  ether  and  carbon 
disulfide  with  frequencies  ranging  from  8.4  to  roughly  2000  cycles 
per  second.  Agreement  between  theory  and  measurement  is  generally 
within  51  with  the  greatest  discrepancies  occurring  near  the  minima. 
Rayleigh  [1890]  and  Mlchte  Smith  [1890]  were  apparently  the  first 
to  use  the  theoretical  relation  as  a  means  of  experimental  deter¬ 
mination  of  T  ,  and  it  has  become  one  of  the  standard  experimental 
procedures.  For  more  recent  developments  and  further  references 
see  Brown  [1936]  and  Tyler  [1941], 

Solutions  for  standing  or  progressive  interfacial  waves,  ana¬ 
logous  to  those  considered  in  section  14  <f  ,  can  be  found  by  appli¬ 
cation  of  the  same  methods.  Since  the  analysis  ts  similar  we  give 
only  the  relation  between  a  and  a-w  .  If  the  two  fluids  fill  the 
whole  space,  with  their  interface  at  -u  *  o  ,  then 


v 

O  = 


r.-/; 

f.  *  ft 


? 


^  r  S 

f 


If  the  lower  fluid  is  of  depth  dp  the  upper  of  depth  d2 
interface  at  •  c  ,  then 


(24.14) 
and  the 


<j  \  f  f'  ~  ft  +  T 

f,  Co iy(x,  V'l  ci,  f-  /vk,  c~ij^ 

In  both  (24.14)  and  (24.15)  er1"  o  if 


(24.15) 


(24,  16) 


thus  the  notion  may  be  stable  even  when  the  lower  fluid  is  less 
dens©  than  the  upper  one.  This  is  not  true  in  the  absence  of  sur¬ 
face  tension,  as  Inspection  of  (14.29)  and  (14.30)  shows. 


The  analogue  of  the  next  example  of  section  14  /  is  somewhat 
more  complex,  for  two  surface  tensions  are  necessary.  Let  T  be 
the  surface  tension  at  the  free  surface  y  *  O  ,  and  T,  z  that  at  the 
interface  ^  *  -  cLf  -  let  the  rigid  bottom  be  at  ^  *•  -/v  -  -  d,  -  . 

Then  the  relation  analogous  to  (14.31)  is 


cr'1'  'lT 

f-vjL  +  (M.l7) 

r"w  fi  +  ~^p  a'*  )  u +  ft( '  *  + 

*  ^  '•  1  >z« 


(24. 17) 
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The  assumption  of  d,»  ao  no  longer  results  in  any  notable  simpli¬ 
fication  of  the  equation.  However,  one  may  show  that  the  solutions 
tr’/j'KV  are  always  real,  and  that  they  are  positive  if 

f,  +  Aw^!  (24.  18) 

This  Is  the  same  condition  for  stability  as  was  found  in  (24.16) 
(and  is  still  necessary  as  well  as  sufficient).  Nuch  of  the  rest 
of  the  pure  gravity-wav©  analysis  of  section  14  cf  may  be  carried 
through.  Thus,  if  <r  is  the  larger  and  the  smaller  root  of 
(24.17)  for  a  given  <v*v  ,  then  one  may  establish  the  inequality 


<  --  <r(/ 1  -  <  Si  ✓ 


(24. 19) 
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If  fj|  and  *7/2  the  profiles  of  the  free  surface  and.  interface, 

respectively,  then  one  find#,  analogously  to  (14,34), 


V 

1  ' 


GO*L  m  dt 


(24.20) 


again,  it  follow  front  (24.18)  that  this  ratio  Is  positive  for  the 
larger  and  negative  for  the  s nut ller  of  the  two  roots  of  {24.17), 
The  discussion  of  the  nature  of  the  motion  associated  with  the 
root  CT^  may  be  taken  directly  from  section  14  S  ;  however,  the 
upper  bound  for  the  velocity  •*£  of  a  progressive  wav©  >f  wav© 
number  4n  is  now  given  by 
t 


nr1 

2k 

•v 

/tH 


±  c  J./I+Z 


T  i 


<  J_  h  -f.  J—  <(vv  ct  i  ^  ^"l(  ^  ^  P  M  /,  _  .  . 
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Let  us  turn  next  to  the  situation  in  which  the  two  fluids  are 
moving  and  look  for  possible  steady  motions.  Assume  each  flu’d 
to  move  to  the  left  with  mean  velocity  C£  and  let  l  *  4  2  , 

be  the  complex,  velocity  potentials.  We  again  look  for  solutions 
in  the  form  [cf.  (14.36)] 

('j  )  m  ~  ci  J  t  ft  (]  )  /  '24 . 22) 

where  is  assumed  small  with  respect  to  6;  J  .  Then  the 

linearized  boundary  conditions  corresponding  to  (14.37)  are 


'l(X)  Ivk  f.  (x  -  -  o;  -  f 

r  P-e  j  fi(x  tv  o)t  C  f J  (jt  +  t  o)  j  -  -L!  Zt  j^f,  (x 


(24.23) 
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If  we  look  for  a  steady  motion  of  the  form 
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then  substitution  in  (24.23)  yields 


and 


(24.24) 


(  f,  C, 7  +  cL\)  «.  (  f,  -  +  r < 24 . 25 > 

The  last  equation  will  not  have  a  real  solution  for  4»vi  assuming 
f,  7  ,  unless 

Vfi-  r*  ,r  £(fic.l+  f^i)1  .  (24.20) 

There  are  then  two  solutions  of  the  form  (24.24).  The  effect  of 
surface  tension  may  be  seen  more  clearly  if  one  graphs  each  side 
of  (24.25)  and  finds  the  intersections,  if  any.  It  will  be  shown 
in  section  26  that  this  type  of  motion  is  unstable  if  I  Ct  "  / 


becomes  too  large. 

Singular  solutions.  The  methods  used  in  section  13  Jf  for 
finding  source-type  solutions  can  generally  be  extended  to  take 
account  of  surface  tension.  Aside  from  the  algebraic  complica¬ 
tions  the  chief  difficulties  are  associated  with  selecting  the 
proper  boundary  conditions  at  infinity.  For  a  stationary  source 
of  pulsating  strength  one  may  still  impose  a  radiation  condition 
as  in  (13.9)  and  obtain  the  correct  solution.  However,  for  the 
steadily  moving  source  of  constant  strength  the  proper  choice 
is  no  longer  clear.  Although  it  is  possible  to  fall  back  upon 
arguments  based  upon  considerations  of  group  velocity,  they  are 

'fill 
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aot  completely  convincing,  so  that  it  see**  .safer  to  formulate 
first  an  initial-value  problem  which  can  yield  either  of  the  two 
cases  mentioned  above  as  a  limit  when  t  -*•*>,  First  me  give  the 
velocity  potential  for  a  source  of  variable  strength  t  le  , 

moving  on  an  arbitrary  path  (4. c{tj)  The  potential -function 
^  must  satisfy  the  same  conditions  given  on  p.  90  except  that 
2)  is  now  replaced  by 


fttfaoi}'*)*  *  f  ?•„■)(*, 0,},t)-O.  (24.27) 

There  is  no  special  difficulty  Involved  in  finding  fj>  .  For 
infinite  depth,  it  is  as  follows: 


+) 


TiD 


(i  )  ^ 

>■  Or) 


(24,28) 
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yet -4  dr  ^.(r/ur-llt J‘i}  Je  '  TJ9(4cfi(f)) 


where  we  have  written  7  for  T/f  .  One  ..  y  similarly  find  the 
function  analogous  to  (13.53)  by  replacing  ^  by  <*)■&  + 7  6  ,  Know¬ 
ledge  of  these  functions  allows  one  now  to  repeat,,  at  least  in 
part,  the  considerations  of  sections  22  c *  and  22  fb  . 

For  ,a  stationary  source  at  (&■/■&,')  with  strength  An  0 *>crt  , 
the  velocity  potential  may  be  easily  derived  from  (24.28),  It 
is  as  follows. 


HA>yt  hf)“H  +  Tt  +  2<3“1' 


7 — - — —  e  4 1  ^  a>0 ert 
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(24.29) 


/  0"  yi 


•f 1  Z  /l  4KH.. 


. . . . I  lllllli!i 


12* 


where  ito  is  the  real  solution  of  a-  ’v  ^4  f  ^ /t  .  It  th©  fluid  is 
of  depth  h,  then 


„  T4  +  j  4  -  o~l’o‘>44v  l*.  L 


1  t"  (t.  kr 

1*.  +  CT 


(24.30) 
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+  ZJTm  T4>-^3  |  t  ■ 

*  J%  £  ■%*  y  ^  3t  ^ 

1  *j>  jr 

where  4.*>  is  the  real  root  of  T  4  +  4  —  cr  co^k.  44  —  O  • 

The  velocity  potential  for  a  source  moving  in  the  direction 
OX  with  constant  velocity  4/t©  may  also  be  obtained  from  (24.38) 
by  a  suitable  limiting  procedure,  although  the  computation  is  some¬ 
what  more  tedious.  In  a  coordinate  system  moving  with  the  source 
it  is  as  follows  for  4  *  «o  : 

Cf(X>  j)  *  ?  "  “T  (24.31) 

n 

—  t?  ^ 

+  U-}  p /a*  ^t0 

JT  l  j  <}+T  4t 
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^  ,  ,, 
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where 


Arc  coi.(4^TV  /ito  h/ 

^  o  ;/  4^r'  *  ut 


and 


•j  4t.4  C<*>*  -  4  ^ 
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On©  may  easily  show  that 

2>  uc'tP'K  &,(&)  4  /  |  /7* '  4  kx.  (^)i  Ax.  c*S‘J"/ T  ' . 

As  T'-*"  0  it  i«  than  evident  that  the  integral  involving 

vanishts  and  that  (24.31)  reduces  to  (13,36). 

On©  nay  carry  out  an  asymptotic  investigation  of  (24,31),  or 

of  9j(  ,  along  the  lines  of  (13.3t>)  and  following.  However,  the 

aaalyeie  is  considerably  more  complicated.  The  behavior  of  the 

as  lonows.  ror 
■i 


wave 


pattern  is  roughly  as  follows.  For  U,f  /4<j  T  1 


for  all  oC 


and  the  disturbance  is 


h)  is  0(H  ) 

chiefly  local.  There  is  a  constant  C  >  1  such  that  when  1  <  ' 

X,  C  the  wave  pattern  is  a  superposition  of  two  set®  of  waves 
corresponding  to  the  two  roots  4,  and  4g_.  Those  corresponding; 
to  *4. |  are  capillary  waves  which  precede  the  source  and  bend 

I  fl 

around  it  so  that  their  crests  eventually  make  an  angle  - 
with  the  X  -axis.  The  gravity  waves  corresponding  to  &■/  behave 
similarly  except  that  they  follow  the  source  and  are  longer.  If 
<1* / 4 <%T'  7  C  ,  a  second  angle,  say  ,  appears,  where  ^  d 

There  are  now  three  sets  of  waves.  Those  associated  with  4j behave 
as  described  above.  The  gravity  waves,  however,  consist  of  both 
transverse  waves  spanning  the  angle  between  ±  (j  Xt  '$tJ  and  diverging 
waves  which  now  lie  in  the  wedge  bounded  by  -  J  +  ^  and  T  f  iA» 
and  its  reflection.  One  will  find  a  sketch,  in  Lamb’s  Hydrodynamics 
[1932,  p.  470]  which  was  computed  for  the  similar  problem  of  a 
moving  pressure-  point,  the  so-called  "fish line  problem".  A  precise 
value  for  the  constant  C  does  not  seen,  to  be  known.  The  free  sui- 
face  t|  may  be  computed  from 
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In  spit©  of  the  general  complexity  of  the  asymptotic  analysis  of 
(34.31),  *t  is  relatively  easy  to  find  the  asymptotic  form,  of  ^ 
directly  ahead (  X  -0°)  and  directly  behind  (of  -180°) : 

e>  0  :  (24.32) 
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here 
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r  f 

and  we  assume  7  4T  ^  .  One  may  see  rather  clearly  the  effect 

upon  4,  and  4-.  r  of  varying  7~  and  *  by  finding  them  as  the 
intersection  of  the  graphs  of  and  >a„l  ^  . 

There  is  no  special  difficulty  in  finding  source  solutions 
for  two-dlmens i ona 1  motion,  and  the  asymptotic  behavior  is  of 
course  easier  to  determine.  The  related  problem  of  a  moving  con¬ 
centrated  pressure  is  treated  in  Lamb  [  1932,  ££  270,  l  ] ,  For  this 
problem  a  paper  by  DePrfma  and  Wu  [1957]  is  particularly  instruc¬ 
tive,  for  they  obtain  the  solution  by  first  formulating  the  initial, 
value  problem  and  then  finding  the  limit  as  t  •©  .  in  addition, 
they  analyze  the  form  of  the  surface  for  large  but  finite  values 
of  t  , 
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2ft •  laves | In  i  viscous  flu, Ad, 

If  one  abandon h  the  assumption  of  a  perfect  fluid  with  irro- 
tetlonal  mat  ion,  on#  loses  at  the  saw#  time  many  convenient  and 
powerful  mathematical  tool--  from  potential  theory  and  the  theory 
of  functions  of  a  complex  variable.  However,  the  simplifications 
introduced  by  the  infinitesimal -wave  approximation  are  sufficient 
to  allow  obtaining  a  number  of  solutions  of  interest,  most  of 
which  have  been  known  for  many  years.  However,  discovery  of  errors 
In  early  work  has  resulted  In  several  recent  papers.  Furthermore, 
in  connection  with  the  theory  of  stability  of  interfaces  the  sub¬ 
ject  has  again  attracted  attention,  this  work  will  be  summarised 

I  In  section  26,  One  will  find  general  expositions  of  many  of  the 

fundamental  results  in  Lamb  [  1932,  $  $ 348-351  ] ,  and  i . evich  ('1952, 

pp.  467-497],  Longuet-Higgias  [ 1953b]  gives  a  valuable  discussion 
of  the  perturbation  procedure  and  carries  through  certain  second- 
order  computations. 

Subject  to  the  limitations  of  the  approximation  one  can  find 
solutions  for  periodic  standing  waves  in  fluid  of  both  infinite 
and  finite  depth  with  a  free  surface,  at  the  interface  of  two 
different  fluids  in  which  either  may  have  a  fixed  horizontal  plane 
as  its  other  boundary,  and  at  the  interface  and  free  surface  when 
two  different  fluids  are  superposed,  the  upper  one  having  a  free 
surface.  In  all  cases  the  presence  of  surface  tension  may  be 
admitted,  By  making  use  of  such  solutions  together  with  Fourier 
analysis  one  can  find  the  solution  to  the  Cauchy -Poisson  initial- 
value  problem  (cf .  Sre tenskil  [1941]). 


In  general,  in  the  investigation  of  standing  waves  one  is 
particularly  Interested  In  two  things,  the  effect  of  viscosity 
upon  the  relation  between  wave-length  and  frequency,  and  the  rate 
of  decay  of  amplitude.  As  an  alternative  to  examining  the  rate 
of  decay,  one  may  instead  assume  that  a  space-and  time-periodic 
pressure  has  been  applied  to  the  free  surface  and  determine  the 
rate  of  transfer  of  energy  necessary  to  maintain  a  steady  oscilla¬ 
tion. 

One  may  still,  as  for  perfect  fluids,  combine  standing-wav© 
solutions  which  are  out  of  phase  in  order  to  form  progressive  waves 
In  a  coordinate  system  moving  with  the  waves  the  wave  system  will 
be  stationary  but  the  motion  will  not  be  steady  for,  as  a  result 
of  viscosity,  it  will  decay  unless  a  periodic  pressure  distribution 
is  moving  with  the  waves  and  doing  work  upon  the  fluid.  Fourier 
analysis  may  be  used  to  obtain  the  fluid  motion  resulting  from  an 
arbitrary  moving  pressure  distribution.  Indeed,  one  need  not  re¬ 
strict  oneself  to  a  pressure  distribution  but  may  also  Include  a 
distribution  of  shearing  stress  at  the  free  surface.  If  a  pres¬ 
sure  and  shearing  distribution  of  localized  extent  is  moving  over 
the  fluid  the  dissipation  of  wave  energy  in  viscosity  will  show 
up  in  a  diminution  of  amplitude,  as  one  moves  away  from  the  pres¬ 
sure  area,  which  is  more  rapid  than  for  a  perfect  fluid.  Such 
problems  have  been  investigated  by  Sretenskii  [1941,  1957]  and  by 
Wu  and  Messick  [1938].  The  latter  include  the  effect  '-f  surface 
tension  and  make  a  particularly  thorough  study  of  the  behavior  of 
the  solution;  they  restrict  themselves  to  two-dimensional  motion. 
One  should  keep  in  mind  that  if  the  fluid  is  of  finite  depth  it  is 
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no  longer  wquivnleat  to  formulate  a  problem  in  which  tho  pressure 
distribution  1®  fiwMl  and  the  fluid  wove®  with  a  eonstaat  wean 
voloelty. 

Instead  of  attempting  to  construct  a  steady  progvooolvo-wave 
solution  by  mean®  of  a  moving  pressure  distribution,  on#  may  in¬ 
stead  assuato  that  tho  progressive  waves  have  been  somehow  initiated 
and  then  study  their  rat®  of  decay  with  distance  from  the  wave- 
maker.  (This  i®,  of  course,  closely  related  to  finding  the  decay 
with  time  of  an  initially  given  progressive  wave.)  Studies  of 
this  nature  have  been  made  by  Biesel  [1940]  and  Carry  [1066] »  who 
Investigated  especially  the  effect  of  the  bottom,  by  Urawll  [1952], 
who  investigated  the  effect  of  side  walls  for  infinite  depth,  and 
by  Hunt  | 1952],  who  combined  the  two.  Dissipation  with  distance 
when  no  walls  are  present  has  been  treated  by  Dmitriev  [1953]  in 
connection  with  the  theory  of  the  wave -maker.  A  point  of  physical 
interest  in  these  studies  is  the  relative  contribution  to  dissi¬ 
pation  of  shearing  motion  near  the  surface,  near  the  bottom,  near 
the  walls,  and  within  the  fluid.  Case  and  Parkinson  [1957]  have 
studied  the  damping  of  standing  waves  in  a  circular  cylinder  of 
finite  depth,  making  use  of  the  linearized  equations  of  this  sec¬ 
tion;  their  experimental  data  seem  to  confirm  the  theoretical 
predictions  when  the  cylinder  walls  are  sufficiently  polished. 

The  fluid  motion  resulting  from  a  submerged  stationary  source 
of  pulsing  strength  has  been  derived  by  Dmitriev  [1953]  for  two- 
dimensional  motion  and  infinite  depth.  Sretonskii  [1957]  has 
carried  through  the  calculations  for  steady  motion  of  a  source 
in  three  dimensions.  Unfortunately,  the  source  function  is  not 
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now  a«  useful  a,  tool  for  constructing  solutions  to  spsclal  boundary- 
vii lue  problems  a*  it  is  for  perfect  fluids.  In  particular,  one 
can  no  longer  satisfy  the  proper  boundary  condition  on  a  steadily 
Moving  body  by  means  of  distribution*  of  source®  and  sinks,  as  was 
possible  in  section  20  jl  .  On  the  other  hand,  distributions  of 
pulsating  sources  may  still  be  used  to  satisfy  the  linearised  boun¬ 
dary  conditions  on  certain  types  of  stationary  oscillating  bodies. 
Thus,,  if  the  motion  is  such  that  the  linearized  boundary  condition 
specifies  the  velocity  normal  to  the  surface  together  with  zero 
tangential  velocity,  then,  a  source  distribution  may  prove  useful. 

For  example,  the  wave-maker  problems  formulated  in  (19.26)  and 
(19,31)  may  be  treated  in  this  fashion;  Dmitriev  [1953]  has  done 

this. 

A  fundamental  assumption  of  the  preceding  rem  .rks  is  that  the 
motion  is  laminar.  Such  an  assumption  seems  to  be  in  harmony  with 
the  assumption  of  small  motions  which  is  made  in  deriving  the  equa¬ 
tions  of  the  present  section.  However,  the  possible  occurrence  of 
turbulent  motion  in  progressive  waves  has  been  reported  by  Dmitriev 
and  Bonohkovskaya  [1955]  who  found  experimental  evidence  for  it  near 
the  surface,  where  the  vorticity  wps  highest.  The  photographs  in 
Figure  8  do  not  seem  to  show  any  evidence  of  it,  but  this  may  re¬ 
sult  from  special  circumstances  of  the  experiments  Bowden  [1950] 
has  essayed  a  theory  baaed  on  von  Harman's  similarity  hypothesis; 
further  references  are  given  there.  In  the  case  of  steady  free- 
surface  flow  in  a  channel  the  importance  of  turbulence  in  modifying 
the  mean -velocity  profile  is  almost  obvious.  However,  investiga¬ 
tions  have  been  confined  to  the  necessary  modifications  of  the 
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shallow-water  approximation  and  will  be  discussed  elsewhere. 

Finally,  we  note  that  much  of  the  theory  given  below  for  a 
constant  surface  tension  T  can,  in  fact,  be  extended  to  a  more 
general  surface  condition.  This  is  Indicated  in  Lamb  [1932,  £ 351] 
and  carried  out  by  Dorrestein  [ 1951)  in  some  detail  for  infinite 
depth.  He  includes  compressibility  of  the  surface  film,  hysteresis 
and  a  "surface  viscosity". 

25  of  .  Linearized  equations  and  simple  solutions;. 

<* 

The  linearized  equations  and  boundary  conditions  have  alreadj 

been  derived  in  section  10.  For  a  stratified  fluid  with  interface 

at  Aj  m  c  the  zeroth -order  equations  are  given  in  (10.2),  the  first- 

order  in  (10.3).  For  a  single  fluid  with  free  surface  they  aro 

given  in  (10.4).  It  Is  customary  and  convenient  to  combine  the 

zeroth-  and  first-order  equations.  Thus,  if  in  (10.4)  we  let 
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p  -  yj  Y  £  pO)  and  v  ■  t  v  ,  then  the  equations  become 
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(23. 1) 
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One  may  clearly  combine  (10.2)  and  (10.3)  In  the  same  way.  In  order 
to  obtain  the  proper  equations  in  a  coordinate  system  moving  to  the 

A 

right  with  velocity  At,  ,  one  need,  only  replace  by  -At*  /6jc  . 


ib1 


TfllliT . ?| . !| 


The  standard  procedure  for  solving  the  equations  is  to  repre¬ 
sent  the  motion  as  a  potential  flow  plus  a  rotational  flow  and  to 
determine  the  pressure  from,  the  potential  part.  Thus ,  let 


where 


and  let 


satisfy 
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(25.2) 
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is  from 

the 

second  equation  in  (25.1)  that 

l/^'must 

y{rl)  -v  £ 

(*■) 

i  V 

(25.5) 

The  relation  between  V  (r>)  and  v_(fl)  is  established  through  the 


boundary  conditions.  In  the  several  examples  treated  below  the 
motion  is  two-dimensional.  However,  there  is  no  difficulty  in 
principle  and  not  much  additional  algebraic  complexity  in  solving 
the  analogous  three-dimensional  problems.  The  essential  simplifi¬ 
cation  in  two  dimensions  is  that  the  components  of  V  ^  may  be 
expressed,  as  a  consequence  of  the  continuity  equation,  in  terms 
of  a  single  function  V  : 

(n.) 
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It  then  follows  easily  from  (25.5)  that 
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Standing  wave;:  - 1  iif tnite  depth.  Wo  shall  try  tc  find  a  solu- 
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t  ton  to  the  equations  which  has  a  profile  of  the  form 

(4*1  je  +M  j  .  (25.8) 

If  such  a  solution  exists,  the  nature  of  A  (i)  will,  of  course,  be 
of  ©e p* *c  ai  Interest.  We  take  ^  and  V  of  the  forms 

&  m  F  /  ir)  t«-3  ( 4u  x  i ^  /  Y  «  G  ■k*'-'  ( /l**'  ■*■*■•<}'  (25.9) 


Equation  (25.7)  then  implies  that 
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where 
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4  1  A  *A> 

'W,  +  — 


(25,  11) 

Neither  JL  nor  u>  need  be  real.  The  form  of  ^  Is  further  deter¬ 
mined  by  Atf>  -  O  and  its  relation  to  Y  through  the  third  boun¬ 
dary  condition  in  (25.1).  It  must  be 


uot 


4*  ~  (a-e  1  +■  to  £  ^7  )  C  C^D  C  ma.  +  <=< 


) 


( 25. 12) 


If,  as  usual,  we  require  the  motion  to  remain  bounded  as  ^ -0 o, 
we  must  take  4  =■  o  .  If  4  has  a  non-vanishing  real  part,  which 
we  assume  for  the  present,  we  may  without  loss  of  generality  take 
it  to  be  positive.  Hence  one  must  have  A  -  o  .  It  follows  from 
the  third  condition  of  (25.1)  that 

pi,  t- 

(K*  G  - - -  *  (25.13) 

<£■  'Vvv 

Substitution  in  the  formula  for  and  integration  with  respect 
to  t  yield 
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♦7  .  c  7" - e  £4*  f*ujt  +  *)-A.e  (25.14) 

f  Z.  V-vtt. 

Finally,  on©  must  substitute  into  th©  dynamical  boundary  condition 
In  (25.1).  Thor©  p»  is  computed  from  (25.4)  with  ^«o  .  For 

future  use  w©  retain  th©  external  pressure  distribution  , 


which  we  take  In  th©  form 


f®  »  j°.  '  U 9»  (/WcJ l  -hoc  )  , 


(25.15) 


where  p,  may  be  complex.  The  boundary  condition  yields  an  equa¬ 
tion  relating  /  and  *w  : 


V  (t  Wj  4  V  "w’-'t  +  *] 


t  »>  /  .  r1  l  P*  2  h"** 

4.  f  <J  **v  f-  /  4*t  *  -  ♦«  . — r-~  * 


f  <- 
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or,  by  making  use  of  (25.11),  an  equation  relating  a*  and  my  : 


+  ?_  ,v,lv) 1  -  4  vL: t  ^  +■  + t  W5 - . 

Consider  f  irst  equation  (25.16)  with  “  <2  a©d  let 

« ,  i  i/  _  +-r^i 

J  ^  1  ^  --- 

> 

Then  (25.16)  takes  the  dimensionless  form 


(25. 17) 


(25. 18) 


(  ^  1  +  I )  v -  4  5  +&•<?. 

An  examination  of  this  equation  shows  that  two  of  its  roots  are 
always  complex  with  negative  real  parts.  These  roots  are  dis¬ 
carded  since  the  corresponding  motion  would  not  die  out  as 
in  fact,  we  explicitly  assumed  earlier  that  /  has  a  positive  real 
part.  (Mote  that  if  we  had  made  th©  other  possible  assumption, 
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i>«»*  that  4  had  a  negative  real  part,  the  resulting  equation 
corresponding  to  (25.16)  would  have  had  roots  with  positive  real. 


part,  again  to  be  discarded.)  The  other  two  roots  have  positive 
real  part.  Whether  or  not  there  is  an  imaginary  part  depends  upon 
the  value  of  K  .  There  is  a  critical  value  *  .581  such  that 

if  ^  M c  the  two  allowable  solutions  are  both  real.  If  K  >  %c, 
the  solutions  are  complex  conjugates.  Let  the  two  complex  root* 
of  positive  real  part  be  denoted  by  ^i±i  ,  Then  one  may 
establish  that  ^  683.  When  the  two  admissible  roots  are  real, 

both  of  them  lie  between  0  and  . 

One  may  find  the  values  of  ^  associated  with  the  two  admissi¬ 
ble  roots  from  (25. 11) .  If  they  are  both  real  <  ^c),  then 
u> •  -v (v*1'- JLl)l  0,  In  this  case  the  motion  is  critically  damped  and 
the  initial  configuration  of  the  surface  gradually  subsides.  This 
occurs  for  a  given  vn  if  V  is  large  enough.  On.  the  other  hand, 
no  matter  how  small  -v  is,  it  also  occurs  when  Avu  is  large  enough, 
l.e.,  for  very  small  wave  length.  If  the  two  admissible  roots  are 
complex  (K  7  tec),  then 
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and 


(25. 19) 


One  may  establish  that  f  -4,  Av»l  v  f\  /+*  J  .  534 ,  8o  that  this  motion 

consists  of  damped  standing-wave  oscillations.  The  larger  /k.  i®, 
the  more  quickly  it  is  damped. 
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Because  of  tbs  relative  complexity  of  equations  (25.16)  and 
(25.17),  it  is  convenient  and  leads  to  more  perspicuous  results 
to1  find  the  relation  between  CO  and  4K.  in  the  two  limiting  cases 
oT  small  and  large  viscosity.  First  consider  the  case  of  small 
viscosity.  If  in  (25.17)  one  lets  V  •*  O  ,  one  regains  the  rela- 
t  ton  o>  '-j'tK.-T.v  of  (24.9);  let  «%  -«p*H.4T WS.  However,  if  one 
re  tains  all  terns  of  the  first  power  in  o'  ,  (25.17)  becomes 

(a/' 4  4  v  +  (25.20) 

which  has  roots 


-  I 


Hence  the  surface  profile  can  be  described  by 


7»  A.i 


4  t  )c9c(m*.x  +■*). 


(25.21) 


(25.22) 


To  this  order  of  approximation,  the  frequency  er,  is  related  to 
44%.  as  In  a  perfect  fluid,  but  the  amplitude  is  gradually  damped. 
To  have  some  idea  of  the  orders  of  magnitude  involved  In  the  damp 
ing,  one  should  consult  the  table  on  page  402  where  the  row  ir0 
gives  computations  relevant  to  this. 

In  order  to  find  the  behavior  for  large  v  ,  divide  equation 
(25,17)  by  and  expand  the  term  [l+w/^ ]  in  a  series. 

If  one  i  a  tains  only  terms  in  and  "V  ,  the  resulting  equa¬ 

tion  leads  to 


i  to14  4  tmU'v  to  4  2(4  <*k,  4-  T  )  *  O , 

I  he  two  solutions,  bofli  real  and  negative,  ate  approximately 
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Her®  1**1  <  and  hone*  is  the  more  important  root  Inasmuch 

ns  it  represents  a  slower  damping  of  the  mot loo.  As  Is  pointed 
out  by  Lamb  [ 1932,  p,  628],  the  root  corresponds  to  a  value 
of  {,  only  slightly  less  than  /Vw  ,  so  that  the  motion  is  nearly 
irrotat tonal.  It  should  also  be  noted  that  by  different  method* 
of  analyzing:  (25.17)  for  large  v  one  may  obtain  somewhat  different; 
coefficients  for  u^. 

In  the  preceding  analysis  it  was  assumed  explicitly  that  £■ 
had  a  non-vanishing  real  part.  If  Z  is  pure  imaginary,  t  *  i  t  , 
another  family  of  solutions  exists.  It  is  now  convenient  to  write 

<jj  and  Y"  In  the  forms 

/  *")  njt  ,  . 

P  *  a  C-  '  €■  iOo  (  Ma  X  t  at  /  J 

Y-  U  CfioZ'y  1-  d  (25.25) 

where 
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The  motion  Is  thus  a  purely  subsiding  one.  The  boundary  conditions 
determine  the  following  relations  between  G,  ,  C  ,  and  <L  : 

a  •  C  ~~ — —  ,  d  «  C  ■—  .  [•»*(*£•  4*$+  +  T ^1 

(25,27) 

All  real  values  of  Z.  are  now  admissible.  The  surface  profile  is 
given  by 
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(25. 2d) 


The  two  sets  of  solutions  may  now  be  used  to  investigate  the  develop 
itent  of  an  initial  disturbance  (cf.  Sretenskii  [1041]). 
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Forced  standing  waves.  *e  may  apply  equations  (25. 1G)  or 
(25.17)  to  answer  the  following  question.  Suppose  that  d*v  is 
given.  Can  we  determine  f>*  in  such  a  way  that  a  steady  standing 
wave 

-  L  <T  t 

*J  *  A0r  c-*>  ( jc  *  <x )  (25.29) 

<.) !  prescribed  frequency  CT  is  maintained?  From  (25.17)  fo*  is  then 
determined  b; 
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v  e 

If,  for  small  viscosity,  one  discards  terms  higher  than  the 
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If  we  take  cr*'*  +  7'!*v1',  the  frequency  obtained  from  perfect- 

fluid  theory,  the  necessary  pressure  distribution  becomes 


—  ■  r-  f  , 

-  4  cr  /vn.  /4  „  t  £  6<^o  iX  +  pc' 


(25.32) 


Thus  the  pressure  must  lead  the  surface  displacement,  by  a  quarter 
of  a  period. 

Standing  wave3-finlte  depth.  If  the  fluid  is  of  depth  h,  the 
analysis  is  similar  to  that  above,  but  yields  expressions  of  much 
greater  complexity.  The  functions  and  ^  may  be  shown  to  have 
the  forms 
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where  again 
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Let 
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Then  C  and  ^  are  related  by 
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(  25  #35) 


The  relation  between  C  and  ro  corresponding  to  (25. 16)  bee ones 
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and  the  surface  profile  is 


7  “  r- —  (c  L  +■  dl  €  /L  ^  (W  x  ■*-<*  ) , 
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(25.37) 


The  formulas  become  more  perspicious  in  the  case  of  small  vis¬ 
cosity  and  no  external  pressure  and  exhibit  the  Importance  of  the 
presence  of  the  bottom.  If  in  (25.36)  one  sets  p„  »  O  and  retains 
only  terms  of  order  y~  ,  and  -v  ,  the  following  equation  re- 

sul ts ■ 
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On®  itjr  solve  thin  equation  by  expanding  10  in  povara  of  v*  , 

u>  m  UJ»  t  U,)t  I*1  '+■  twj.  u  -f  .  .  .  t 

p. 

substi fating  in  (25.38)  and  keeping  only  tern*  An  v  »  *v  ’  and  v  . 
The  tern  Independent  of  v  yields  to,  mt  ic”,  ,  where  <x„  is  given 

by  (25.21)  and  is  the  frequency  for  an  lavlacid  fluid  To  the  order 


of  accuracy  consistent  with  (25.38),  one  finds 

,,  w  C - -  ,  t  «  *  onL  fon.4  ♦  cmi,  z1hk-  | 

to  »±  t  er»  -  ( /  ±  l  )i  2<x,  >  to*tdL  v  %ssaL™ . - - j . . . . 

*t  ,/Vk/  ft  ”*  1 

(25.39) 


The  firnt  two  terns  were  given  by  Hough  [1897].  The  correct  ex¬ 
pression  (25.39)  was  first  given  by  Biesel  [1049]:  Hough  had  given 
. 2av.' v  for  the  last  tern  but  he  apparently  made  an  error  in  calcu¬ 
lation,  for  (25.39)  was  derived  Independently  of  Blosnl’n  work  and 

* 

has  also  been  checked  by  Carry  [1956]  (Basset's  analysis  [l888i,  p. 
314]  overlook*  the  terns  in  V%) . 

The  formula  (25.39)  should  be  compared  with  (25.21),  the  cor¬ 
responding  formula  for  infinite  depth.  There  the  effect  of  vis¬ 
cosity  f  ,,ts  only  with  the  first  power  of  V  .  The  dissipation  of 
energy  in  the  body  of  the  fluid  is  evidently  of  less  importance  than 
in  the  vicinity  of  the  bottom.  When  two  fluids  are  superposed,  a 
similar  phenomenon  occurs  in  the  neighborhood  of  the  interface  [cf. 
(25.44) ]. 

Standing  waves -stratified  fluids.  Consider  now  the  situation 
in  which  a  fluid  typified  by  ?,  and  /+,  fill®  the  space  /-]  <o  and 
another  typified  by  Tt  <  Sf  and  the  space  ^  y  o  ,  The  equations 
to  be  satisfied  in  the  two  fluids  and  at  their  interface  are  given 
in  (10.3).  The  method  of  solution  is  analogous  to  that  used  for  a 
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single  fluid,  lavmr,  separate  functions  4 t  »  *  »B<1  ft  * 

arc  needed  for  the  lover  and  upper  fluid*.  For  a  stand  lag -wave 


solution  they  may  toe  taken  in  the  fern 


<f,  rn  CL,  *■  ^  cm  (m.  x  tot.)  f  X  -  4$V  €**  '  e  #  ^  d**-  (to**  tot,)/ 


(25.40) 


<fr  *  at  &  *  £  4'  ^  O-o  (-ptgit  f»<)  ^  f-t,  *  4-,  €•  '  e  *  ^  W  ♦*"«*)/ 


•here  we  aisum  both  and  to  have  positive  real  part*.  AS 
t,  ,t2  and  /*n.  «r©  related  by  the  equation 


w  ■  y. 


(t'-.uC)-  *.(C -***)■ 


(25.41) 


Subs  ti  tut  ion  of  (25.40)  In  the  various  boundary  renditions  at  >j*  ^ 
gives  four  homogeneous  equations  relating  CL,  ,  ,  4o,  ,  and  *^a  . 

The  determinant  of  the  coefficients  set  equal,  to  aero  yields  another 
relation  between  Uj,  ,  t/  ,  ^j.  and  dn*. : 


f,  +  fv  )  Uj  vrr.-fjj  tv*  t  T  ^14*1/  J 


t  ^  CO  <Wv  +*  /( /t  ,  ♦•'V  *t/C/  £ ,  j  m  Q 


(25.42) 


In  the  limiting  case  of  small  viscosity,  (25.42)  gives 


i>j 1  +  JL&±- 

vZTtf/TJl 


,,**•  +  (fj~  ,0. 


/>  r. 


(25.43) 


This  has  the  approximate  solutiona,  when  the  coefficient  of  <o'H 


is  small  relative  to  the  last  term, 


t  Mf  t  i 

11>iv!y*ii********  *******  **j*mj^j^]fi|p]~|‘^$*i^^pflj]^pi|Hj|^|***^||j||^  ^ 

fft/*a7(M.44) 
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where  0i  is  the  psrfcci-fluld  frequescy  gives  is  equation  (84.14). 
This  solution  was  first  Hives  by  Harrison  ( 1908] .  The  most  sign! 
f least  physics)  fact  about  (25,44)  when  coepared  with  (25.21)  that, 
to  the  order  oi  approx  last  ion  considered,  the  latter  shows  a  rate 

t 

of  decay  proport  lost  1  to  4*,  v  and  no  influence  of  viscosity  on  the 

frequency,  whereas  (25.44)  shows  a  rate  of  decay  and  an  alteration 

of  the  frequency  proportional  to  (la  a  dleensional  sense). 

The  greater  importance  of  vlecoaity  for  stratified  fluids  nay  'toe 

ascribed  to  the  different  boundary  condition  at  the  interface. 

Harrison  computed  the  wave  velocity  and  Modulus  of  decay  (tine  re- 

-f 

quired  for  the  amplitude  to  decrease  by  a  factor  e  )  for  an  air- 
eater  interface  at  17°C  (  f,  •  1,  ft,m  0.00129,  ^  *  .0109,  ^  *  0.139, 

T  -  74  in  c.g.s.  units).  In  the  following  table  reproduced  Iron 
Harrison *s  paper  V.  ,  %  and  4r  sre  the  wave  velocities  neglecting, 
respectively,  both  surface  tension  and  viscosity,  viscosity,  and 
neither;  t*  ,  T  ,  are  the  moduli  of  decay  taking  account  of  the 
water  viscosity  only,  a  water-air  Interface  without  surface  tension 
and  a  water-air  interface  with  surface  tension. 


lave -length 

(cm  ) 

i 

10 

100 

1000 

14  (cm /sec) 

12.48 

30.46 

124.79 

394.82 

24.90 

40.05 

124.81 

394. 62 

nT 

24.89 

40.04 

124.81 

304.62 

7 - 

*"-»  l(p 

1.162'* 

1*56,2" 

3  hrs  12*39.4" 

321  hrs  5*40" 

r 

1.125" 

1*34.1" 

1  hr  21*40  6" 

36  hrs  50  *3 6" 

— 

L  c 

!. ,  106" 

1*34.0" 

1  hr  21*40.3" 

36  hrs  50' 3 4" 
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A  striking  aspect  la  the  apparent  importance  of  the  air-water  In¬ 
terface  in  damping  loaf  eaeee  and  almost  total  lack  of  influence 
on  wav©  velocity  (the  latter  fact  is  obvlotie  free  (85.44)), 

For  very  large  viscosities  the  roeulta  are  analogoua  to  those 
for  a  single  fluid.  The  two  values  of  w  analogoua  to  those  in 

(25.24)  are 


"T+’>r . 


2**.  y^,+/*v  ft  t  fu 


(25.45) 


The  analysis  of  the  roots  of  (25.42)  for  general  values  of  V, 

and  -v,  is  difficult.  However,  it  has  been  carried  through  by 

§*• 

Chandrasekhar  I  1955,  especially  pp.  170-173]  for  the  special  situa¬ 
tion  v,  -  Uj,  and  T*  0  •  In  this  case  l,  1  /Z  •  Th*  b*h“vlor 

ls  similar  to  that  described  for  (25.17)  except  that  the  critical 
value  Kc  separating  a  steadily  decaying  motion  fro.  an  oscillatory 

decaying  one  is  now  a  function  of  (f>,  +  M  •  Thi®  VlllU#  <actU>lljr 

a  different  one  since  he  chooses  a  different,  parameter)  is  tabulated 
for  a  variety  of  density  combination*.  Further  analysis  of  <25.17) 
may  be  found  in  a  paper  by  Hide  (1955]  and  Tchen  [ 1050b). 

Kusakov  [1944]  has  carried  through  an  analysis  similar  to 
Harrison's  when  the  upper  fluid  is  of  depth  & t  ,  the  lower  of 
depth  L,.  However,  the  results  do  not  seem  to  be  consistent  with 
Harrison’s  (or  those  above)  when  and  k «.  become  large.  This 
problem  has  also  been  considered  by  Hide  [1955],  but  with  an  approxi¬ 
mation  that  neglects  the  viscous  boundary  conditions  on  the  walls. 
Harrison,  in  the  same  paper,  haw  treated  also  the  problem  when  the 
upper  fluid  is  of  finite  depth  and  with  a  free  surface.  We  shall 


ii 


mrnmn . . . .  ituttHtidiiHii*  HhiiHihi**  m  v  ih-i«i  » *■*  h*  i « *«  »i  i«m  ■*-!**  *  ►•♦ohii-hi  m  *  m  » -  *«'i.  *  i  *  mi 


i  mu  II II  III  M  HI  Ml! . INN W1 *1  M II  "I'l  1 1  I  ll'H'l  'I . .  #  mflWpUh Ml SIP* . 


4  014 


not  reproduce  the  results  except  to  rtwrk  that  hi*  ca*putatioas 
•  ho*  *  very  marked  influence  or  tho  damping  of  a  thin  layer  of 
fluid  of  a  lightly  different  donmlty.  Tho  of foe t  of  a  variable 
surface  tension  upon  eave  motion  la  iaveatigated  briefly  in  Lamb 
[1932,  $  351  ]  and  at  **««*  length  in  Levich  [1952,  pp.  477-490). 

Pulaina:  stationary  source,  Daitrlev  [1953]  has  derived  the 
form  of  the  functions  f  and  XY  and  the  surface  profile  in  the  pres¬ 
ence  of  a  submerged  source  of  pulsating  intensity  -<?  cosert,  le 
■hall  give  here  only  his  expression  for  the  surface  profile  and  an 
asymptotic  expression  for  large  distances  from  the  source.  Let 
the  source  be  located  at  ( 0 ,  -  kt)  and  let 


L  «  It.Jff 


J?  1  *-*/?  1  r)-V 


'■  2T 


ST.’p 
tj  I  «y 


The  surface  profile  is  then  represented  by 


’  Ji.  t 

Qpltrt  ^  (  I  ~  Z-t  X 

7  '  &-*—  f  J  iIJlYYiH*  +X  h)h  -  X*jn(X-tj 


-  L  X  —  y  j  v 
£  CiOXAO-A 


(25. 46) 


4  A  K 


Q  f0+  H>°  1  ) 
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26. 


Stability  of  free  surfaces  and  interfaces 


In  this  section  we  wish  to  examine  the  circumstances  under 
which  a  small  disturbance  of  a  free  surface  or  of  an  interface  be¬ 
tween  two  fluids  will  Increase  in  magnitude  with  time.  The  energy 
for  this  increase  may  come  either  from  available  potential  energy, 
e.g.  if  the  lower  fluid  is  lighter  than  the  upper  one,  available 
kinetic  energy  in  the  case  of  flowing  fluids,  from  forced  motion 
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of  solid  boundaries,  or  possibly  sow*  other  source  such  as  a  given 
pressure  distribution  over  a  free  surface.  Surface  tension  and 
viscosity  nay  be  expected  to  have  a  stabilising  effect,  so  that 
special  interest  attaches  to  the  study  of  their  Influence.  We  shall 
use  the  nature  of  the  energy  source  as  a  convenient  on*  for  separat¬ 
ing  classes  of  problems,  even  though  not  every  situation  falls 
clearly  Into  one  of  the*. 

Since  the  boundary  conditions  and  equations  which  w*  shall  use 
for  the  Mathematical  analysis  have  been  linearised,  following  the 
assumption  that  the  disturbances  sre  snail,  one  cannot  expect  the 
predictions  of  the  theory  to  be  valid  quantitatively  much  beyond 
the  initiation  of  an  unstable  motion.  However,  a  great  advantage 
in  the  use  of  linearized  theory  is  that  sn  arbitrary  initial  dis¬ 
turbance  can  be  analyzed  into  Fourier  components  and  the  behavior 
of  individual  components  examined  separately. 

26  cx .  Interface  be tween  stationary  superposed  fluids . 

Following  our  earlier  notation,  let  us  identify  quantities 
referring  to  the  lower  fluid  by  the  subscript  1  and  to  the  upper 
fluid  by  2.  Let  a  sinusoid  distrubance  of  wave  number  4n>  exist 
at  the  interface.  Consider  first  the  esse  of  perfect  fluids  with 
no  surface  tension.  Then,  if  both  fluids  are  infinitely  deep,  the 
relation  (14.28)  must  hold.  If  f,  7  f t  ,  the  standing-wave  solution 
of  section  14 1  obtains.  However,  if  J,  <  it  ,  then  or  <:  o  and 
cr  is  imaginary.  Let  ,  i.e.  <*>  *  — 4~ 

Then  one  must  replace  an(o'^*T)  in  the  of  that  section  by,  say, 
UjJL  Out  •  The  profile  of  the  free  surface  is  then,  according  to 


CIO. fl),  given  by 


7  *  A  a*.  4tv  jc  c^M.  u>  ^ « 

*  (26. 3) 

The  amplitude  of  the  initial  corrugations  of  the  surface  evidently 
increases  very  rapidly  with  time,  and  the  solution  is  a  valid  ap¬ 
proximation  for  only  a  limited  time  interval.  The  nature  of  the 
disturbance  need  not  have  been  restricted  to  Sc*v  m  jl any  function 
$,)  satisfying  A<f -  O  would  have  yielded  the  same  be¬ 
havior.  Equation  (26.1)  atm  holds  if  the  two  fluids  are  bounded 
below  and  above,  respectively*  by  and  )jm  ®*cept  that  to 

is  given  by 


UJ 


ft  -  f, 


Jxi£l 


(26.3) 


C#44*  VM..  ii  i,  i  t  A  |  Jj;  -V  f f 

The  surface  is  still  unstable,  but  the  rate  of  growth  of  the  ampli¬ 


tude  is  slower. 

Effect  of  surface  tension.  Let  us  now  suppose  that  surface 
tension  acts  at  the  interface.  Then  the  relation  between  &  and  *** 
given  in  (24.14)  or  (24.15)  must  hold,  and  a  standing-wave  solution 
is  possible  even  if  f* /  ,  provided  that  (24.16)  holds,  i.e. 


f, 


(26.4) 


Thus  the  interface  is  stable  under  small  disturbances  of  sufficiently 


small  wave  length.  If  the  unequallty  in  (26,4)  is  reversed  and  we 
again  set  then  (26.2)  holds  once  more  and  the  solution 

is  unstable.  However,  the  value  of  is  less  than  that  when  T-  o  , 

so  that  the  rate  of  growth  of  the  disturbance  is  retarded.  It  is 

t, 

also  clear  from  the  form  of  the  relationship  between  u)  and 
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that  there  la  a  eave  number  for  ehlch  $  that  le  the  rate  of 

growth  o2  the  disturbance,  in  a  maximum.  If  both  fluid*  are  of 

infinite  depth  this  mode  of  maximum  instability  occurs  when 


*  (k  -  fJi/s-- 


(26.5) 


The  effect  of  finite  depth  of  the  fluids  is  to  displace  the  posi¬ 
tion  of  the  maximum  to  higher  values  of  (smaller  wave  lengths) 
but  a  precise  calculation  requires  solving  a  transcendental  equa¬ 
tion. 

Effect  of  viscosity.  The  influence  of  viscosity  in  stabilis¬ 
ing  interfacial  disturbances  has  been  the  subject  of  s  number  of 
recent  papers,  in  particular  Bellman  and  Pennington  [1954], 
Chandrasekhar  [1955],  Hidef 1955]  and  Tchon  [1956].  The  relevant 
equation  relating  U>  and  4n.  is  now  (25.42).  Because  of  the  high 
degree  of  this  equation  it  is  not  easy  to  give  a  complete  dis¬ 
cussion  of  its  admissible  roots.  However,  it  is  easy  to  establish 
that  if 

(  f,  '  k  )  $  +  T  aJ:  <  O  ,  (26.7) 

then  (25.42)  has  a  positive  real  root  W*  satisfying 

9*4-r**%  (16.g) 

Thus  the  presence  of  viscosity  does  not  alter  the  conditions  for 
instability,  as  the  presence  of  surface  tension  did,  but  it  does 
have  a  stabilizing  effect  in  that  the  rate  of  growth  of  a  distur¬ 
bance  is  slower. 

In  order  to  show  the  existence  of  a  positive  root  under  con¬ 
dition  (26.7),  one  can  write  (25.42)  in  the  form 


'Ii|i|'lli  1+11 'll . . . . . .  ll  i|i  IMIi  Ill'll  HWIUlWI  iHtlNI  l<|HIH|i  i| 


4*0? 


(  +  [f  i  m  f  i  *"*v  r  ^  *  **  A  W  **V  . — - 


At ti**v  +i«,  4  f  MifHA/AiL 


ip*** 

and  sketch  as  functions  of  -o  the  curvaa  represented  by  the  two 
sides  ot  the  equation  (remembering  that  and  /tare  functions  of 
c o  ),  The  statement  above  then  follows  easily  from  the  fact  that 
both  curves  are  continuous  and  the  one  represented  by  the  right- 
hand  function  starts  at  the  origin  like 

-  2  ♦**(/* *  v*» ) w 


and  goes  to  -  »o  in  the  fourth  quadrant,  behaving  as  uj~->  oo  like 


„  vk  \ZTTTa/u 

UJ  AW  »j— —  ■ — 1 - 

5. A  ♦  v  f» /t. 


/V* ♦  v  V’ 

A  more  elaborate  discussion  of  the  roots  is  given  by  Bellman  and 
Pennington  [1954). 

The  behavior  of  u>„  as  a  function  of  Avt,  in  the  interval  defined 
by  (26.7)  and  in  particular  the  mode  of  maximum  instability  has 
been  investigated  by  the  authors  cited  earlier.  Chandrasekhar  has 
computed  the  curves  to.fW)  for  V,  -  Vz  ,  T-  o  and  a  number  of  values 
of  (ft  ■  ft  )j(?L  t  ?  i  )  •  Hide  has  recomputed  these  by  an  approximate 

method  and  then  applied  the  method  further  to  a  fluid  of  finite 
depth  with  a  continuous  density  variation  f.  .  Tchen  has 

devised  a  different  method  of  approximate  computation  and  Includes 


the  effect  of  surface  tension.  Figure  33,  which  is  chiefly  quail- 

v 


tative,  shows  the  variation  of  to 
interval  of  instability. 


as  a  function  of  *w  in  the 


FIGURE  33 

Accelerated  fluid.  If  the  whole  system  of  fluid  la  being 
accelerated  In  the  y-direction  by  a  constant  amount  -  ?i  »  then 
the  relative  motion  in  a  moving  coordinate  system  is  the  same  as 
if  the  system  were  at  rest  and  <jj  had  been  replaced  by  j  +  fi  ,  as 
is  immediately  evident  from  equation  (2. 15).  With  this  change  *he 
reasoning  of  the  preceding  paragraphs  still  applies.  This  fact  was 
pointed  out  by  G.  I.  Taylor  [1050]  who,  on  the  basis  of  At,  formu¬ 
lated  the  following  rule  (neglecting  the  influence  of  surface  ten¬ 
sion,)  i  If  the  fluids  are  being  accelerated  in  a  direction  from  the 
more  to  the  less  dense  fluid,  the  interface  is  stable;  in  the  con¬ 
verse  case  it  is  unstable.  Experiments  carried  out  by  Lewis  [1950] 
for  large  accelerations,  about  50g,  confirm  Taylor's  observation 
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and  the  predicted  initial  rata  of  growth,  Taylor' a  paper  gave  rise 
to  a  nurr.be r  of  other#  treating  various  aapecta  of  the  Instability 
of  accelerated  interface*.  In  addition  to  thoee  cited  in  the  last 
paragraph,  we  mention  Ingraham  [1954],  Plesset  [  1954 J,  lirk'hoff 
[1954],  and  Lay ter  [1955]  but  shall  not  summarize  the  content®. 

The  effect  of  an  imposed  acceleration  oscillating  in  magnitude  will 
be  discussed  in  section  26  ^  . 

26^6.  Interface  between  moving  fluids. 

Consider  the  situation  in  which  the  fluid  occupying  the  region 
u  <  t)  ftj  >d)  is  moving  to  the  left  with  velocity  -  t  (~  ci)  ,  and  sup¬ 
pose  that  a  small  disturbance  exists  near  the  Interface.  If  we 
suppose  that  the  fluid  is  perfect  and  the  notion  in  each  fluid  ir- 
rotational,  then  we  may  describe  it  by  the  velocity  potentials 

4 i  (*•  y  <  J 1  *  )  “  ”  C  i  x  +  fl(J' 1  i  1 1  )  ■  (26.9) 

4e  shall  assume  c,  f  Ct  . 

The  kinematic  boundary  condition  at  the  interface  may  be 
written,  after  linearization  appropriate  to  tho  assumption  of  a 
small  disturbance,  in  the  form: 

tyUf'I't)*  +  (26  10) 

The  dynamical  boundary  condition  (3.9)  yields  the  following  generali¬ 
zation  of  (10.8): 

(r  ( f*.(  fit  mT(*lxx+  7}))  k  -«• 

J ( 26 .11) 
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If  *i  in  el  iminated  to#1  twin  (36.10)  and  (30.11)#  on*  find* 


(36.13) 


/ 


C  |  —  C  i 


Tlh w ’<3. 


Let  ua  now  restrict  our  attention  to  two-di sene i one 1  notion 
of  fluids  bounded  above  by  y  bt  and  below  by  '4,  .  and  lot 
the  Initial  dioplaeeaeat  be  *j/(x,o).  Then  fron  (15.2)  we  know  that 
the  subsequent  notion  nay  be  resolved  Into  harmonic  progressive 
waves  moving  to  the  right  and  left.  It  will  be  sufficient  for  our 
purpose  to  examine  a  single  component  of  the  spectrum.  Hence, 
we  look  for  a  solution  in  the  form 

All-  X  ~  ^ 

f'  -  a,  /Hv(^  ivuj  4  / 

(  (  a*,  x  -  <J~i  ) 

<j-\  *  ut  c&*L  /wv.  i  a-j~  *«■*)  e 
It  follows  from  (26.10)  that  (c,  -  Cv)  7*-- 


(26. 13) 


Hence 


h  s  -H-i —  fa,,  i-u»L  ^  b{  +  <3-t  tuX  -Vi*  4t,  1  C  '  ^  ^ 

l  C,~CVL 


It  then  follows  from  (36.10)  that 


(26. 14) 
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U«X  M*  4  ,  +  <*-»,  4V.  4- 1  Q,  I  4Ka  ^  ^ 


C-,  - 


-£a 


C-4a^*» 


Tt* 

(27.15) 


Substitution  of  (26.13)  in  (26.12)  and  use  of  (26.15)  yields  the 
following  relation  between  c r  and  : 
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(26.16) 


The  solution  may  bo  expressed  *•  folio** : 

.21  .  f-i  $>  -2  *  f'1  ^  Ct*^  ^  4i  y  + 

**  ?,  4wu4v  ~ 

(/,-/,)  J  * r«K  ~^7~  fthult-L* t^L,(8»a9) 

.f/  C#^4v  ^  ^ ^  ^ »_  fillip,  4v.  ^  ^  (f{  Ca44a  4*"*  ^  4^  C.|^) 

It  is  evident  from  the  form  of  the  term  under  tbo  radical  that  cr 
cannot  be  real  unless 


■f  Taw  >  (c, 


■<») 


$<_fi  0*44*, » n*  ‘A-/  c*44*  4»a 

f,  cWt,  4%*  I  l. 


(26. 17) 


It  Is  thus  evident  that  there  are  no  real  solutions  unless  the  left- 
hand  aide  is  positive  and  that  there  may  even  then  exist  an  inter¬ 
val  of  wave  numbers  for  which  the  disturbances  is  unstable  (If  both 
o  and  r  are  aero,  such  a  velocity  discontinuity  is  always  un¬ 
stable).  If  one  assumes  f,>  fz  ,  the  minimum  value  of  the  left- 
hand  aide  :1s 


2  | }(f,-f<)fT 

and  occurs  for  a**?  ~(ft  ~  ft)  J j]~  .  Since 


(26.18) 


l 


f,  f’,  C-&44,  iv.  C&4C,  m>-  k  f ,  fi.  (26.19) 
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the  disturbance  will  be  unstable  for  some  wave  numbers  whenever 
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One  may  conclude  from  (26.16)  that  the  horlsontal  walla  hare  a 
destabilizing  effect  In  the  sense  that  wave  numbers  which  are 
stable  for  Infinitely  deep  fluids  may  become  unstable  modes  In  the 
presence  of  walla.  For  an  air-water  interface  the  right  side  of 
(26.20)  Is  about  (646  cm/sec)2.  The  corresponding  wave  length  is 
1.71  cm:  If  the  water  is  at  rest  (c, -  °)  ,  then  the  wave  velocity 

is  0,84  cm/sec  in  the  direction  of  the  wind . 

Let  us  suppose  that  Ct  and  ar®  keth  positive,  i.©»  that 
both  fluids  really  do  flow  to  the  left.  Then  it  follows  from 
(26.10)  that,  if  the  roots  are  real,  one  of  the*  is  always  negative 
and  thus,  fro*  (26.13),  represents  a  wave  moving  along  the  interface 
in  the  direction  of  the  stream.  The  other  will  propagate  upstream 


if 


(f, 


f. )  -i  +  r  /*w  >  f,  c, 
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4k  i  J  ^  Cr%,  /•'W  >4.  ^  . 


(26.21) 


otherwise  also  downstream. 


An  investigation  along  the  above  lines  of  the  stability  of  an 

interface  between  flowing  fluids  was  first  given  by  Kelvin  { 18V 1 ] . 

•! 

Similar  treatments  with  additional  information  may  be  found  in  many 
texts,  especially  Lamb  l 1932,  *$232,  268]  and  Rayleigh's  Thepryof 
Sound  [Cambridge,  1926,4366].  Kelvin's  intention  was  to  try  to 
predict  the  minimum  wind  velocity  which  will  cause  a  small  distur 
bance  on  smooth  water  to  increase  in  amplitude,  and  to  find  the 
unstable  wave  lengths.  The  predicted  minimum  velocity,  roughly 
650  cm/sec,  is  much  higher  than  the  obmerved  minimum,  which  is  about 
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100  cii/see.  An  evident  objection  to  the  analysis'  above  is  that 
viscosity  of  both  air  and  eater  has  been  neglected,  Since  this 
alter#  In  an  essential  way  the  behavior  of  the  fluid#  near  the 
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interface.  It  Is  not  surprising  that  the  prediction  Is  not  accurate. 
One  should  not  espoct  confirmation  except  In  circumstances  In  which 
it  la  possible  to  show  that  the  effect  of  viscosity  in  confined  to 
n  neighborhood  of  the  interface  small  with  respect  to  the  niitliaun 


wave  lengths  considered.  The  subject  of  wind  generation  of  waves 
is  still  in  an  unsettled  state.  One  will  find  surnames  of  the 
present  status  in  the  article  by  H.  U.  Boll  in  vol.  48  of  this 
encyclopedia,  especially  pages  703-717,  and  also  In  a  critical 
exposition  by  Orwell  {1958].  A  summary  of  some  of  the  work  In 
fti#  USSR  on  wave  feneration  is  included  in  Shulelkln  [1956], 

The  inclusion  of  viscosity  in  the  analysis  above  leads  to  a 
somewhat  more  difficult  development  than  In  the  case  of  standing 
waves.  An  exposition  of  the  present  achievements  In  this  theory 
will  be  omitted;  they  consist  chiefly  of  papers  by  Vuest  [1949] 


and  Lock  [1981,  1954). 


2 6  |f .  Ver t 1 call.y  oscillated  basins. 


Let  S  denote  the  wetted  surface  of  a  basin  and  F  the  water 
surface  when  the  basin  is  at  rest.  We  shall  suppose  that  the 

n 

basin  is  being  oscillated  in  the  y-dlrection  according  to  some 
given  law,  which  may  be  specified  by  giving  v9(l)  ,  the  velocity 
of  a  point  of  the  basin.  It  will  be  most  convenient  to  describe 
the  motion  of  the  fluid  in  coordinates  fixed  in  the  basin,  these 
will  be  denoted  by  *,  y,  ^  .  We  shall  assume  the  oscillations  and 

a*  -1,  .  tv-  •  — —  •••■  .  . . . * I 


I"MhMi|<i|I' 


f 


4j?t 


and  the  resulting  notion 


to  bo  of  sMll  nap!  Itudo  so  that  wo  any 


I  linearise  tbo  equations  and  boundary  conditions. 

If  <p  la  the  velocity  potential  for  the  notion  relative  to 
the  basin  and  the  profile  of  the  surface,  both  In  coordinates 
fixed  in  the  basin,  then  it  follows  easily  from  (*.17)  that  the 
1  only  necessary  change  is  to  replace  <j  by  f  +  tr.  in  the  boundary 
conditions  at  the  free  surface.  They  become: 


<26.22) 

+  ft  (■*'<>,),*)  «  7VV)  ,  T'-  T/f  (26.23) 

On  the  basin  walls  one  must  have 

-  0  -  S  (26. 24) 

We  wish,  as  usual,  to  investigate  the  character  of  the  motion  of 


the  fluid. 

The  problem  formulated  above  is  elsarly  related  to  the  problem 
considered  in  section  23^  .  However,  the  resulting  motions  are 
quite  different.  Rayleigh  [1883]  appears  to  have  made  the  first 
theoretical  investigation  of  this  problem.  More  recently  It  has 
been  studied  by  Moiseev  [1053,  1954],  Benjamin  and  Ursell  [1054], 
Schultz-Grunow  [1955]  and  Bolotin  [1956].  Moiseev's  analysis  is 
the  most  general  in  that  the  only  restriction  upon  the  basin  shape 
is  that  it  should  allow  construction  of  a  Green's  function  for  the 
Neumann  problem;  surface  tension  is  not  taken  into  account. 

Benjamin  and  Ursell  restrict  themselves  to  basins  in  the  form  of 
a  vertical  cylinder  with  horizontal  bottom,  but  include  the  effect 


4 1<> 

of  surface  tension.  However,  at  the  intersection  with  the  wails 
they  assume  a  90°  angle  of  contact  with  the  free  surface.  This  is 
In  contradiction  with  the  observed  behavior  of  fluids  but  simplifies 
the  mathematical  treatment.  In  spite  of  this  shortcoming  it  seems 
desirable  to  include  the  effect  of  surface  tension,  and  this  will 
bo  done  below.  Bolotin's  paper  considers  a  modification  for  vis¬ 
cous  damping.  The  treatment  below  follows  closely  that  of  Benjamin 
and  Ursell. 

Let  the  basin  be  of  depth  k  ,  let  C  denote  the  intersection 
of  the  walls  with  the  plane  *o  ,  and  let  ft  be  a  normal  to  the 
wall  at  a  point  of  C  .  Then,  from  (26.22)  and  (20.24)  it  follows 
that  f  *0  ,  or  -  const,  at  each  point  of  C  ;  w©  take 

this  constant  to  be  zero,  thus  assuming  a  90°  contact  angle  with 
the  wall.  It  then  follows  from  (26.23)  that  ^ .  <?  . 

Let  ^  (4,  j )  be  a  set  of  functions  harmonic  in  the  region 
bounded  by  the  basin  and  the  plane  and  satisfying  (26.24), 

and  such  that  (*>  ot])  form  a  complete  set  of  orthonormal  functions 
in  the  area  of  the  (X,})  -plane  bounded  by  C  .  Then  <^Lx,o,  }  ,i )  , 

'V  ( x  ’  aa<1  *  7  J  5.  can  each  *>e  expanded  in  series  in 

(A  (x>  °'  ))  •  Th®  ®  *Pan®ion  of  $  (x,  O'},*)  determines  <j>  (■*,  )/}>+) 
ms  a  series  in  )  •  In  the  case  at  hand,  when  the 

basin  is  a  vertical  cylinder,  one  may  separate  variables  as  in 
section  12  ck  and  construct  a  set  ^  in  the  form 

jZ  t )  ~  (x < P 

\  ,  (26.25) 
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ITh©  ©  a v ft  L iji^hs  4n^  are  determined  by  the  boundary  condition  on  the 
contour  C  ,  namely 

Let  the  expansion  for  be  written  in  the  form 


lU.}.t>-ZaL(V  ?t^}) 


(20,27) 


Then,  by  differentiating  (26.27)  and  using  (26.26)  one  gets 


(26.28) 


then 
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and,  from  (26.22) 
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Hence 
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( 26 . 29) 


Now  substitute  (26. 27) -(26.29)  in  the  remaining  boundary  condition 
(26.23)  : 

• Z  [( y  +  v".  j  A.  t  7”  /Vw  ^  +  -  <3.  ^  c&tA,  K-. 

Since  the  ! £  are  orthogonal,  we  may  set  each  coefficient  of  (f>^ 
equal  to  zero.  With  the  special  choice 


tC  *  C  Ced  <T  t 


(26.30) 
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the  following  set  of  differential  equations  determine  the 


+[( +  1  **,4  )  -^4  4  +C4k1 ^  WX  coo  crtj 0*^(4 )  m  c 

(26.31) 

1 f  we  net 

T  '  t  rrt  '  (?*v''4t  T **4.  )  *  ^K.  > 

(26.32) 


where  OJj^  is  the  frequency  of  free  oscillations  in  the  mode  <wv^ 
when  the  basin  is  fixed,  then  (26.31)  takes  one  of  the  standard 
forms  for  the  Mmthieu  equation: 


Xr1 


a,t  + 

1* 


(26.33) 


Of  particular  interest  in  the  present  context  is  the  behavior 
of  the  solutions  ^  as  't t  or  if",  becomes  large.  It  is  known  from 
the  theory  of  differential  equations  with  periodic  coefficients 
that  a  pair  of  fundamental  solutions  can  be  given  in  the  form 

e  Q  (?)  t  z  Q(~z)/  (26.34) 

where  Q  is  of  period  IT  ,  unless  is  an  integer.  In  the  latter 
case  there  exists  a  periodic  solution,  of  period  JT  if  y>’  is  even 
and  of  period  2  J  if  odd,  and  another  independent  nonperiodic  solu¬ 
tion.  The  coefficient  j*.  will  be  a  function  of  the  parameters  , 
£4  and  it  is  particularly  pertinent  to  the  present  investigation 


to  know  for  what  regions  in  the  fag)  -plane  ju.  has  a  nonzero  real 
part.  These  regions  have  been  investigated  for  other  purposes  and 
■ay  be  found,  for  exaaple,  in  M.  f.  McLachlan's  Theory  and  applica¬ 
tion  of  Mathleu  functions  [Oxford,  1047,  pp.  40,41].  In  Figure  34, 
reproduced  from  Benjamin  and  Ursell,  the  shaded  regions  represent 
the  unstable  regions  of  the  (ft  f  )  -plane  where  T  has  a  nonzero  real 
part.  In  the  unshaded  regions  ju.  is  pure  imaginary  (but  not  an  in¬ 
teger)  and  the  two  solutions  (26.34)  are  bounded  for  all  T  ,  The 


J 


boundaries  between  regions  correspond  to  the  periodic  solutions 
occurring  when  yu.  is  an  integer.  In  the  unstable  regions  the 
periodicity  behavior  of  the  solutions  is  of  two  types.  In  the 


II 
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•ecoftd,-  fourth)...  reg  Ions  is  real  and  the  solutions  (86.34)  are 
functions  of  period  W  multiplied  by  exponent in Is.  Ira  the  first, 
third  .  reg i on*^ •ju, r 4 reel,  rand  the  solution!®  (26.34)  now  be¬ 
come  functions  of  period  2  3T  multiplied  by  exponent  Inis.  In  tern* 
of  t  the  two  set*  of  region*  cor r era pond,  respectively,  to  frequen¬ 
cies  c r  and.  |  O’  . 

For  a  given  mode  of  oscillation  one  must  compute  and ^ ^ 
and  plot(%|^on  the  stability  chart  in  order  to  find  out  whether 
the  inode  is  stable  or  not.  It  seems  likely,  and,  in  fact,  has 
been  proved  by  Moiseev  [1954,  p.  44],  that  for  any  given  values  of 
or  and  c  some  of  the  possible  modes  will  be  unstable.  However, 
the  analysis  above  has  neglected  the  damping  effect  of  viscosity 
and  it  nay  be  supposed  that  the  only  unstable  modes  which  actually 
occur  are  those  associated  with  the  smaller  values  of^fc.  In  any 
case,  as  has  been  emphasised  earlier,  the  analysis  in  only  suitable 
for  describing  the  initial  stages  of  the  motion. 

If  the  frequency  of  oscillation  er  is  equal,  or  nearly  so,  to 

one  of  the  frequencies  cr^  for  free  oscillation  of  the  fluid,  or 

% 

to  a  subharmonic  of  ,  1 . e .  cr  *  cr^  ,  then.  f\ml  ,  or  1a.  , 

and  It  is  evident  from  Figure  34  that  '£■&.)  will  be  in  an  unstable 

region.  If  crrcr^  wlU  lie  in  the  lowest  r®*lon  <md  ®ta“d“ 

irag  waves  with  half  the  frequency  of  the  basin  will  be  generated. 

If  cr  ml  (3^  wAl1  lle  ln  the  second  region  and  the  generated 

standing  waves  will  have  the  same  frequency  as  the  basin.  It  is 
pointed  out  by  Benjamin  and  Ursell  that  an  apparent  discrepancy 
between  experimental  observations  of  Faraday  and  Rayleigh  and  of 
Matthiessen  can  be  explained  by  the  above  remarks. 


1 1  iiimiwwiNf mmm  «n  i  0  Wfff 


4m 

Benjamin  and  Ursel 1  made  an  experimental  Investigation  with  a 
circular  cylinder  in  order  to  determine  by  experiment  the  boundaries 
of  the  lowest  region  of  instability.  The  measurements  provide  i* 
surprisingly  good  confirmation  within  certain  limitations. 

27.  Higher-order  theory  of  Infinitesimal  waves. 

It  is  implicit  in  the  theory  of  infinitesimal  waves  developed 
in  the  preceding  sections  of  this  chapter  that  the  approximation 
given  by  first-order  theory  to  the  solution  of  a  particular  problem, 
assuming  that  one  exists,  can  be  improved  by  including  further  terms 
in  the  perturbation  aeries.  The  solution  of  the  resulting  bound  try - 
value  problems,  at  least  in  the  simplest  cases,  can  be  carried  through 
in  a  manner  similar  to  that  of  the  first-order  theory,  although  the 
computations  become  more  and  more  tedious  the  higher  the  order  of 
approximation.  Nevertheless,  in  view  of  the  interest  of  the  re¬ 
sults,  the  computations  have  been  carried  through  by  a  number  of 
persons  and  by  a  variety  of  methods. 

Stokes  [1849]  was  apparently  the  first  to  make  the  calculation 
for  progressive  waves;  in  fact,  the  method  used  below  in  section 
27  is  not  essentially  different  from  Stokes’  first  method.  Later, 
in  connection  with  the  publication  of  his  collected  papers,  Stokes 
[1880]  added  a  supplement  describing  a  different  procedure.  Rayleigh 
turned  to  the  problem  several  times  [1876,  1911,  1915,  1917]  and 
introduced  still  another  method  of  approximation.  It  should  be 
noted,  however,  that  both  Stokes*  second  method  and  Rayleigh's 
method  are  limited  to  two-dimenaional  irrotational  progressive  waves. 
Rayleigh  [1915]  seems  to  be  the  first  to  have  given  an  adequate 
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treatment  of  the  higher-order  theory  of  standing  waves.  In  addition 
to  then®  classical  pa*\.rs  there  have  been  many  others  extending  or 
laiproving  the  earlier  theory;  sore  of  these  will  be  noted  below. 

In  all  such  computations,  and  indeed  in  the  numerous  first- 
order  computations  carried  out  in  the  earlier  sections  of  this 
chapter,  there  is  the  tacit  assumption  that  there  exists  an  ’’exact 
solution"  which  Is  being  approximated  and  which  can  be  approached 
more  amd  more  closely  by  pursuing  the  selected  method  of  approxi¬ 
mation.  Unfortunately,  it  is  seldom  that  one  Is  able  to  prove  the 
existence  of  an  exact  solution  or  of  convergence  of  the  method  of 
approximations,  and,  in  fact,  Burnside  [1916]  cast  doubt  upon  the 
usefulness  of  the  Stokes-Ray leigh  type  of  approximation  of  periodic 
progressive  waves  of  permanent  type.  Burnside ' s  objection  was 
later  met  by  Nekrasov’s  [1921,  1922,  1951J,  Levi-Civita ’s  [1925] 
and  Struik's  [1926]  proofs  of  the  existence  of  such  waves  for  both 
infinite  and  finite  depth.  However,  the  existence  of  a  standing 
wave  satisfying  the  exact  boundary  conditions  has  not  been  demon¬ 
strated  as  yet.  The  same  is  true  of  the  more  complicated  problems 
considered  in  earlier  sections.  However,  this  mathematical  t  nort- 
comlng  is  possibly  of  no  more  importance  than  the  neglect  in  many 
problems  of  relevant  physical  parameters  such  as  viscosity. 

One  should  bear  in  mind  that  the  higher-order  infinitesimal 

<fi 

waves  considered  below  are  not  the  only  higher-order  approximations. 

The  solitary  and  cnoidal  waves  of  the  next  chapter  bear  a  similar 
relation  to  the  first-order  shallow-water  theory.  In  addition,  in 
the  last  chapter  another  method  of  approximating  exact  waves,  due 
to  Havelock  [1918},  will  be  described. 


In  the  following  we  shell  be  seeking  a  wave  which  •loveii  with¬ 
out  change  of  form,  i.e.  a  progressive  wave  in  the  sense  of  section 
7.  Hence  we  shall  expect  to  be  able  to  represent  ^  and  In  the 
form 

0  -  f(x  -  C  t  ,  y,})  ,  «  .  *j(x  -  ct  )  ,  (27.  1) 

where  C  is  the  velocity  of  the  wave.  It  will  be  convenient  to  re¬ 
present  the  motion  in  a  motion  In  a  moving  coordinate  system,  say 
X-x-ct-  However,  we  shall  henceforth  drop  the  bar  over  the  ^  , 
The  boundary  conditions  at  the  free  surface  are  then  the  follow¬ 
ing  : 

%(*,))  %  })-  +  7.  ^  -  c  ^  ■  *  '  (27.2) 


c  f)1  Yr;'+0-°, 


(27.3) 


where  K,  +■  is  given  by  (3.5’)  and,  as  usual,  .  Surface 

tension  is  being  takeD  into  account  both  for  the  intrinsic  interest 
ol  the  results  and  because  of  an  interesting  phenomenon  which  occurs 


in  the  higher-order  approximations.  We  shall  suppose  that  the  wave 
length  A*  of  the  wave  system  has  been  given,  so  that  c  is 


still  an  unknown  of  the  problem. 

bet  us  now,  as  in  section  10  <x  ,  assume  that  ?  ,  *f  and  C 
may  all  be  expanded  in  a  perturbation  series  in  some  parameter  £  : 

+  e  7  *'  <a7-4) 
C  ~  <T.  t  t  C,  +  Ct  +  •  •  •  . 


iHwimii'i 


4w 


After  substituting  in  (27.2)  and  (27.3)  and  collecting  terns  In 
the  manner  of  section  10  ,  one  obtains  the  following  boundary  con¬ 


ditions  which  must  be  satisfied  successively  by  cf^\  £m 

7  (iJ*  c2  : 
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(27.5) 
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where  all  conditions  are  to  be  satisfied  on  the  plane  aj  ~  o  •  It 
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Is  possible,  of  course,  to  carry  the  approximations  further,  but 
three  steps  are  ample  to  illustrate  the  procedures.  The  solution 
will  be  carried  through  in  outline  through  the  third  order  for  in¬ 
finite  depth  and  through  the  second  order  for  finite  depth.  As  an 
expansion  parameter  we  may  take  t  ^Am,  ,  where  A  is  a  length  de¬ 
termining  the  amplitude  of  the  waves.  The  motion  will  be  restricted 
to  be  two-dimensional. 
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Infinite)  depth.  The  solutions  of  (27. 5)  are  already  known 
from  (13,5).  We  take  the*  in  the  following  for* 


,c//»  C. 


■  C  iOV  /VH  X  ,  *J  *  ^  0*0  JC  t 
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*  ^  0*0  Afth.  X  ,  t*  4*1#  «  #■  4**.  f 


(37.8) 


After  substitution  in  (27.9),  one  find* 

foQ-)  <t» 

4  <T.  *}  C,  Oua,  tn.  j<  —  C.  Lck,  £■♦**<  JC, 
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rtt ) 

Elimination  of  ♦£  yields 
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(27. 10) 
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aa  the  boundary  condition  to  be  satisfied  by  r  .  If  c'  f*0' 
one  cannot  find  a  periodic  potential  function  satisfying  (27. 10) 


Bonce  we  set 


C,  -  o, 


(27. 11) 


A  solution  of  Laplace’s  equation  satisfying  (27.10)  with  ct*o  and 


vanishing  as 


.oc  is  easily  found  to  be 


A}  s  c.  V  r' 
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(27.  12) 
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providing  4u.  #  j/tr1  •  The  corresponding  expression  for  ^  is 


.ft,)  i  <  5  t/w,1  T 1  .. 
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(27.13) 


One  could,  of  course,  add  terms  of  the  form  given  in  (27.8)  but 
with  arbitrary  multipliers.  However,  such  solution®  are  discarded 
since  we  wish  to  allow  only  first-order  terms  of  this  form. 


4  "2.4 


Two  striking  facts  show  up  in  (27.12)  and  (27,13): 

First,  it  surface-  tension  Is  neglected,  <pl?)  vanishes  and  ‘f 
gives  the  velocity  potential  correctly  to  at  least  the  second  order. 

/pit)  <tj 

The  second  fact  is  the  zero  in  the  denominator  in  both  t  and  7  } 
which  shows  that  <^?  and  become  unbounded  as  ^  approaches 
f  $ /C7  '  .  Cif)  may  argue,  of  course,  that  this  simply  shows  that 

validity  of  the  perturbation  method  is  limited  to  smaller  and  smaller 

- - - p 

values  of  A  **  the  closer  one  comes  to  ,  However,  it  seems 

. . . « 

also  to  be  an  Indication  that  near  -  \$/ir'  the  mode  represented 

by  9  u‘>  is  of  the  same  order  of  magnitude  as  that  represented  by 

,o  0 ) 

T  .  That  this  is  indeed  the  case  is  clear  from  an  examination 

of  the  'equation  determining  '  '  and  ^  ^  when  4m,  «  if/tT'  ,  In 

rj>  (l) 

fact,  T  was  not  determined  by  (27.  10)  for  this  value  of  and, 
furthermore,  (27.8)  does  not  give  the  complete  solution  of  (27.5). 

The  solution  with  which  we  must  start  in  this  case  is 


Cf(>mSl  (  w.jt  +CL  t  ^  X  +  4  &  yc&o2*«jtJ 

-mv  L  , 


/  (27. 14) 


where  a  and  -C  are  as  yet  undetermined  constants.  Thus  these  two 
modes  of  motion  are  of  the  same  order  for  M^m  £  7 Tr71.  one  may 
now  substitute  (27.14)  and  the  corresponding  -7  ^  into  (27.9).  By 
reasoning  similar  to  that  used  earlier  in  setting  C,  *  o  ,  we  now 


find 


CL  *  i 


1 
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1  £  “  0  i 


(27. IS) 


There  are  thus  two  possible  first-order  modes  depending  upon  the 
sign  of  Cl  ,  is  now  a  sum  of  terms  with  modes  tnd 

f <u  4 au .g  »  but  will  not  be  given  here.  The  wave  profile,  including 


■edits  through  on  2  a*  jt  «  **y  b«  writ  ton  •*  follow*: 
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(27. 17) 


The  two  signs  in  the  second  solution  coriv  . r  id  rough lj  to  the 
change  of  sign,  occurring  In  the  first  when  k  panne**  through 

. .  Comparison  of  the  two  cures  also  gives  an,  indication 

of  the  limitations  upon  A  +k  necessary  in  the  first  solution, 
namely, 


fA^I  <  h 

1  1  I  $  +  h*,1  r  / 

A  reversal  of  curvature  at  the  center  of  the  wave  trough  for 

v*  <  i  o  /zt'  ►  or  of  the  C1'est  for  ijTzr  ,  will  occur  when 
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The  existence  of  the  singularity  in  the  expressions  for 
and  f^.aa  ftrat  noticed  by  H.rri.on  (1900].  Hilton  [1915] 
mined  the  matter  more  carefully,  found  the  solution!!  (27.17)  and, 
in  fact,  carried  all  approximations  further.  Soji®  of  Wilton's  com¬ 


puted  profiles  are  shown  m  Figure  35.  Although  Wilton  casts 
doubt  upon  the  existence  of  the  solution  (27.17)  with  ♦$,  such 
profiles  seen  to  have  been  observed  by  Kamesvara  Rav  [1920],  How 
ever,  the  matter  apparently  still  awaits  a  thorough  experimental 
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Investigation.  <jk>  also  similar  higher  nodes  mentioned  below. 


Let  us  now  turn  to  the  next  order,  assuming  **v  f  f  ffzT*  .  Sub¬ 
stitution  of  (27.8)  and  (27. 11) -(27. 13)  Into  (27.7)  and  elimination 

(t  ) 

of  yields  the  following  boundary  condition  to  be  satisfied  by 

f(i)  on  y  =  0  : 
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Again  In  order  to  avoid  an  unbounded  solution  we  must  set  the  coe 
fficient  of  w  "o  x  equal  to  zero.  This  yields  a  value  for  cz  : 
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One  may  now  find  a  potential  function  satisfying  (27.20)  and  vanish- 

(p  j  _  * , 

Ing  as  ^  — *>  -  oo  .  The  solutions  for  r  and  ^  are  as  follows: 


f  j'i  _  j?  £o  ***  r'  £?  +  z*ul  r ') _  ^  z 
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(27. 22) 
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for  4t\.  j  & fhT'  .  From  (27.22)  one  sees  again  that  p 

would  vanish  if  surface  tension  were  neglected.  Although  we  shall 
not  carry  through  the  computation,  this  does  not  happen  for 


It  IB  kLbo  evident  thnt  another  singularity  has  appeared  at 


•n  i 


-  i%/lT'  '  In  fact,  when  one  examines  the  reason  for  the  appeer 
ance  of  the  singularities,  it  i«  evident  that  a  node  of  the  for* 

03  will  always  show  a  singularity  at  ei  *  I  f  /***  7*  .  In  each 

such  case  the  reason  is  the  same  as  in  the  situation  discussed 
earlier  with  *  2  :  for  4*v  -  the  proper  first-order  solu 

tion  is  of  the  for*; 
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with  Q,^  to  be  determined  subsequently  (according  to  Wilton  only 
civ  is  not  unique),  Thus  (27.8)  should  be  qualified  by  4u.  t  ^/mJ  • 
One  should  note  that,  although  44m ^  is  getting  small  (and  hence 
large)  as  increases,  the  wave  number  of  th©  second  ft  nut -order 
mode  is  fa,  ^  /T‘  ■  Hence,  on  the  basis  of  the  results  in  section 
25,  one  will  expect  this  mode  to  be  quicVly  damped  for  large  values 
of  'Vv  .  However,  one  may  presume  the  first  few  to  be  observable. 

We  remark  that  these  special  associated  pairs  of  first-order  waves 
always  straddle  the  wave  number  for  minimum  <-0  ,  namely  A*i  . 

The  wave  profile,  velocity  potential  and  wave  velocity  will 
be  given  by 
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c  *  C.  +  A  'tv*  C,  +  A  •'W  ct  +  •  *  ■  .  (27.24) 

To  the  third  order  the  profile  for  pure  gravity  waves  ( T  '*  o  )  is 
represented  by  the  following  function: 
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A  •  A  i*  +  }  A  J  ;  the  velocity  becomes 
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The  velocity  potential  to  the  third  order  is 
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(27.27) 


If  one  sets  a  -  O  ,  then  the  wave  profile  for  pure  capillary  waves 


becomes 
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and  the  velocity 
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For  pur©  gravity  waves  the  approximations  were  carried  to  the  fifth 
order  by  Stokes,  Rayleigh  [1917]  and  others. 

It  is  of  interest  to  compare  the  profiles  represented  in  (27.25) 
and  (27.28).  The  effect  of  including  higher-order  terms  in  pure 
gravity  waves  is  to  sharpen  and  raise  the  crests  and  to  broaden  and 
raise  the  troughs.  For  p;  %  capillary  waves  the  effect  is  just  the 
reverse.  For  combined  gravity-capillary  waves  the  increasing  im- 
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port  a  nee  of  the  second-order  term  near  4*+  *  I  ff'zr  ‘  first  shot 

up  as  a  reversal  of  curvature  at  the  Middle  of  the  flattened  part 


of  the  wave ;  formula  (27.19)  gives  the  condition  for  the  first 
occurrence.  In  Figure  35  are  shown  a  pure  gravity  wave  as  computed 
by  Wilton  [1914]  for  *  .66  (here  A  is  the  amplitude) ,  and  five 
ty-capi 1 lary  waves,  the  last  two  corresponding  to  the  solutions 
7)  j  also  computed  toy  Vi  1  ton  [1919]  for  a  value  of  l‘/f$  m  .  075. 


grav 
( 27 . 


It  should  be  remarked  that  the  value  of  .86  is  much  larger 

than  any  for  which  it  is  possible  to  prove  convergence  of  the 
perturbation  series  and  is,  in  fact,  very  close  to  the  value  of 
/l,1lVlii  for  the  highest  possible  irrotational  wave  of  permanent  type 

(see  st-  T  ion  33i>  •  ,  ran:«iy  0.891. 


Finite1  depth.  When  a  solid  bottom  is  present  at  aj 
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only  necessary  modification  of  the  preceding  analysis  is  substltu- 

ro(t- )  .  i  (// i  t 

t  Ion  of  the  boundary  condition  J  -k,/*P  for  t  o  as  y*~wo 
This  inc i eases  the  computational  labor  by  a  substantial  amount, 


but  otherwise  introduces  no  difficulties.  However,  we  call  atten¬ 
tion  to  the  remarks  on  the  definition  of  wave  velocity  in  section 


the  veloc i ty 


below  is  the  one  defined  there  also  a 


The  wave  profile,  velocity  potential  and  wave  velocity,  in¬ 
cluding  the  effect  of  surface  tension,  are  as  follows,  to  the 

second  order 
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The  velocity  is  the  same  as  la  the  first-order  theory;  this  occurred 
also  for  infinite  depth.  In  contrast  to  the  case  of  infinite  depth, 
the  term  '  7  does  not  vanish  when  1  m  °  .  The  singularity  in  the 
coefficient  of  ee i  Z,4h, x  still  persists  provided  that  4? 

The  earlier  discussion  of  this  phenomenon  is  still  relevant,  and 
a  detailed  one  will  be  omitted  here.  However,  even  if  surface 
tension  is  neglected  in  (27.30),  the  second-order  term  may  still 
become  large  for  small  values  of  •w. 4  ,  as  has  been  empha.s lined  by 
Miche  [1944].  If  one  again  takes  as  an  indication  of  increasing 
predominance  of  the  second-order  term  a  reversal  of  curvature  at 
the  bottom  of  the  trough,  one  finds  that  this  occurs  for 

A  i  Ux~M.  ***  ^  7 

M'h-  /  j  /  7— ;  r . /  (27.33) 
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or  approximately 

,<4  y  -i  ix>viA^c  b-’  'Vvt  hs 

as  given  by  Miche.  The  occurrence  of  this  secondary  crest  when 
Avv  is  small  has  frequently  been  observed.  It  has  been  Investi¬ 
gated  experimental  1 ly  by  Morison  and  Crooke  [1953]  and  by  Horikawa 
and  Hegel  [1059]. 

The  wave  profile  and  velocity  computations  were  carried  by 
Stokes  to  the  third  order  for  pure  gravity  waves  in  fluid  of  finite 
depth.  The  following  expressions  are  taken  from  a  report  by 
Skje Ibreia  [1959] 
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Skjelbreia  has  provided  comprehensive  tables  allowing  easy  compu¬ 


tation  of 


t  '  f  and  many  other  quantities  of  interest,  all  to 
the  third  order. 

Particle  orbits.  A  particularly  Interesting  phenomenon  occurs 
when  higher-order  approximations  are  used  In  the  computation  of  the 
paths  of  individual  particles.  The  equations  which  the  coordinates 
of  a  particle  must  satisfy  are 
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(27.35) 


di  -r~  V 

Since  'P  depends  upon  the  parameter  £  ,  the  solutlonB  X  and  J 
also  will.  We  assume  then  that  X  and  h  may  be  expanded  into 


series  of  the  form 


X(i)-  x..t*,(i)+  ,  te).  i.+  ei.t)*...,  (27  36) 

substitute  them  into  (27.35)  together  with  the  appropriate  expan¬ 
sion  of  in  powers  of  t  ,  and  then  equate  the  several  powers  of 

!r  separately.  This  results  ir?  a  sequence  of  equations  of  which 
the  first  two  are  as  follows 
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The  first  set,  (27.37),  was  already  solved  In  (14.17)  and  (14.18) 
and  to  the  first  order  of  approximation  gave  circular  or  elliptical 


orbits. 

surface 

vanish. 


The  solution  for  higher  orders  is  facilitated  by  neglecting 
tension  and  assuming  L  m  00  ,  for  then  ¥'  and  f  both 
From  (27.8)  one  finds  easily  the  orbit  to  the  second  order: 
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The  circular  orbits  of  first-order  theory  are  now  modified  by  a 
general  drift  in  the  direction  of  wave  motion.  The  total  amount 
of  fluid  transported  per  unit  time  (and  width)  is  z  A  am.  C0  .  As 
the  formula  shows,  this  additional  flow  is  concentrated  chiefly 


near  the  surface. 

When  the  depth  is  finite,  or  when  surface  tension  is  taken 
into  account,  the  orbits  become  more  complicated,  bet 
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The  particle  orbits,  accurate  to  the  second  order,  are  as  follows: 


iJ  3.  (9. 
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The  mass-transport  term  in  x(ti  is  still  present,  and  In  fact, 
persists  to  the  very  bottom.  The  elliptical  orbits  of  the  first- 
order  theory  are  now  modified  not  only  by  the  forward  drift  at  all 
levels,  but  also  by  another  superposed  cyclic  motion  of  twice  the 
frequency.  The  effect  of  this  is  to  make  the  orbits  approximately 
epitrochoidal  (neglecting  for  a  moment  the  drift)  with  a  small 
hump  at  the  bottom  which  in  extreme  cases  can  become  a  cusp  or  a 
loop.  This  behavior  has,  in  fact,  been  observed  by  Morison  and 
Crooke  [1953].  For  capillary  waves  the  situation  is  reversed  and 
a  dimple  appears  at  the  top, 

lhe  existence  of  mass  transport  will  be  reconsidered  in  the 
last  chapter,  where  it  will  be  demonstrated  that  it  is  a  general 
consequence  of  Ir rotational  motion  when  the  exact  boundary  condi¬ 
tions  are  satisfied.  The  theoretically  predicted  monoton ically 
decreasing  forward  drift  with  increasing  depth  is  not  confirmed 
experimental  lv  for  small  values  of  44  ,  say  4  4  <  2.  .  Instead, 
with  respect  to  a  coordinate  system  moving  with  the  mean  velocity 
of  the  fluid,  there  is  an  observed  forward  flow  near  the  bottom 
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and  top  and  a  backward  flow  in  the  middle  portions,  It  in  not 
surprizing  that  the  perfect-f luid  nodal  doas  not  giva  a  good  pra- 
dlction  for  snail  4  -4/  ,  for  the  high  shear  rate  near  the  bottom 
Indicates  that  viscosity  should  aot  be  neglected.  Longuet-Bigglns 
(1953b]  has,  in  fact,  devoted  a  long  nonograph  to  development  of 
the  higher-order  theory  of  waves  in  a  viscous  fluid  and  finds  theo- 

-f 

retical  drift  curves  agreeing  qualitatively  with  observed  ones. 

We  shall  not  carry  through  the  details  here  and  refer  to  Longuet- 
Higglns'  paper. 

Wave  energy.  One  of  the  striking  facts  about  progressive 
first-order  pure  gravity  waves  is  that  the  kinetic  and  potential 
energy  per  wave  length  are  equal  (see  section  15/*>),  This  equal 
division  of  energy  no  longer  holds  when  higher-order  terms  are 
taken  into  account.  It  is  particularly  easy  to  show  this  for  k=<*>, 
for  then  we  may  use  (27.25)  and  (27.27).  The  average  potential 
energy  in  a  wave  length  is 


(27.42) 


The  average  kinetic  energy  is 
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Composite  waves.  Previously  in  this  section  we  have  been  dis¬ 


cussing  a  wave  of  permanent  type  whose  prototype  is  the  first-order 
progressive  wave  of  the  form  « »  A  ci  ^(x-d).  It  Is  natural  to 
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Blesel  [1952]  has  derived  formulas  for  a  composite  wave  with  a 
finite  number  of  components  and  for  k  finite.  He  computes  a  number 
of  quantities  of  interest.  However,  the  formulas  are  very  long  and 


will  not  be 


suiting'  motion. 
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attention  of  several  persons*  Among  these  art  iretemskll  [1911]# 
Bergman  and  thappelear  f  1957) ,  Daubert  [l®§7,  IfM)  in  a  series  of 


notes,  Joins  [IMS)  and  Nomad  is  [  1957 


U' 


Bretenekli  [1959 


$ 


has  investigated,  the  higher -order  theory  of  save  motion  rums l ting 
from  a  moving  pressure  distribution  and  waves  in  a  circular  carnal. 
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waves . 


As  will  bo  evident  bo  low,  the  formulation  of  a  higher-order 


theory  of  standing  waves  is  somewhat  clumsier  than  that  for  pro* 
gross ive  waves  of  permanent  type*  Fart  of  the  difficulty  stems 
from  the  fact  that  on#  necessarily  must  deal  with  one  more  variable 
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The  type  of  motion  we  are  seeking  will  he  represented 


by  a  profile  ^(x,  tj  and  a  velocity  potential  periodic  In 


both  X 
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If  -wo  fix  the  wave  length  X  ”  ^/fc,  >  then  the  period  7  y4r  will 
have  to  be  determined  as  one  of  the  unknowns  of  the  problem.  In 


addition,  we  wish  to  have  the  first -order  standing  wave 
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ti 

of  section  14  oi  serve  as  a  prototype  and  first-order  solution  of 
the  no  re  general  problem*  As  a  further  condition*  w#  shall  suppose 
the  motion  to  be  symmetric  with  respect  to  a  vertical  line  through 
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Rayleigh  [1915]  was  apparently  the  first  to  consider  this 
problem.  It  was  later  attacked  in  an  entirely  different  way,  using 


agr aitg  i  an  c oor d  1  nates ,  by  Sell#1  r» h - 2e nk o v  I  c  b  [  1947' ,  19 9 1  a ,  b ,  19  92  J 


'who  treated  tooth  two-  and  three -dimensional  waves  for  infinite  depth 
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two-dlMnsiOHtl  waves  for  finite  depth,  and  compos i to  eaves  for  In 


finlto  depth.  Peaoey  and  Price  [19M],  following  approximately 


Raylo Igh'o  method,  Carr lod  the  approaimatlon  for  two-dimens tonal 
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otlM  and  1..*^  to  the  fifth  order*  mud  to  tho  second  order  for 
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pans  ions  for  an  infinitely  deep  fluid. 
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the  situation  is 


similar  to  that'  discussed  in.  section  JIT  *C  following  (IT.  13)* 
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Aw  bu  been  antitaH  earlier,  the  approximation  has  beea  carried 


to  tho  fifth  order  by  Penooy  aad  Price 


cover,  it  is  not 


oeceaaary  to  carry  the  approximation  eo  far  ia  order  to  aeo  hm 


Important  features  of  the  notion,  namely  the  sharpening  of  the 


crests  and  flattening  of  the  troughs,  the  absence  of  any  nodal 


points  and  the  decrease  of  frequency  with  amplitude 


One  interest¬ 


ing  feature  doe a  require  carrying  tho  approximation  to  at  least  the 


fourth  order :  this  in  the  absence  of  any  time  during  n  period  when 


the  surface  is  cm 


flat.  In  connection  with  am  experimental 


tost  of  a  predicted  standing  mate'  of  greatest  amplitude -length 


ratio  by  Penney  and  Price,  ffi.  t*  Taylor  [1® 


has  also  provide! 


an  experimental  verification  of  the  correctness;  of  the  theory  In 


an  extreme  case 
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n.  The  method  of  computation  of  or  to  its  including  higher 
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order  terms  Is  the  same  as  that  outlined  nt  the  end  of  section  27 o( 
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The  effect  of  the  last  term  In 


in  easily  seen  to  be  a  nun 11 


wiggle  superposed  on  tho  t tret-order  nt might- linn  trajectories 


discussed  in  section  14  tx  t  except  directly  beneath  the  crests 


where  the  trajectory  ia  nil 11  vertical  but  with  the 'midpoint  some 


what  above  the  eqjulllbrltin  position 
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We  note  that  here  also,  as  in  the  case  of  progressive  waves,  the 


terms  in  and  y 


4vv  L 


0 


21 


Waves  in  a  viscous  fluid 


The  equations  (10. 2) -(10. 4) ,  used  **»  section  26  in  developing 


first-order  theory  of  waves 


a  viscous  fluid,  may  be  con 


i 

aide red  as  the*  first  in  a.  sequence  for  the  determination  of  higher 


order  approximations.  Although 


formulation  of  the  equations 


ppears  to  he  straight  forward,  if  laborious,  the  higher-order 
heory  does  not  seem  to  have  attracted  many  investigators. 
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Harrison  ( 1909 ]  made  a  second-order  investigation  of  progressive 

waves  and  Longue t -Higgins  [1953]  has  recently  nade  an  elaborate 
study  of  both  progressive  and  standing  waves  in  an  attempt  to  ex- 
plain  certain  observed  features  of  mass  transport  velocities.  We 


shall  not  at temp 


summarize 


#r»  the  follow 


Aug  results*  taken  from  Harrison,  may  be  of  interest.  For  the  wave 
profile  to  the  second  order  he  gives  the  following  expression  when 
1  flj  s  mS1 1 1  ( c  f «  (  >2  5 « 21  >55 ) )  « 
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where  CT# " ■«  j4^,  The  effect  of  viscosity,  besides  damping,  is  to 

make  the  leading  side  of  the  crest  steeper  than  the  trailing  side. 
According  to  Harrison  the  average  horizontal  velocity  of  a  particle 


again  for  smal  X  T>  ,  is 
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This 


formula  should  be  compared  with  A  ^  computed  from  (27.39), 


to  which  it  reduces  when  *x)  •  o  . 
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